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Abstract 

Modifications of general relativity provide an alternative explanation to dark 
energy for the observed acceleration of the universe. We review recent develop- 
ments in modified gravity theories, focusing on higher dimensional approaches 
and chamelcon//(i?) theories. We classify these models in terms of the screen- 
ing mechanisms that enable such theories to approach general relativity on small 
scales (and thus satisfy solar system constraints). We describe general features 
of the modified Friedman equation in such theories. 

The second half of this review describes experimental tests of gravity in 
light of the new theoretical approaches. We summarize the high precision tests 
of gravity on laboratory and solar system scales. We describe in some detail tests 
on astrophysical scales ranging from ~ kpc (galaxy scales) to ~ Gpc (large-scale 
structure). These tests rely on the growth and inter-relationship of perturba- 
tions in the metric potentials, density and velocity fields which can be measured 
using gravitational lensing, galaxy cluster abundances, galaxy clustering and the 
Integrated Sachs- Wolfe effect. A robust way to interpret observations is by con- 
straining effective parameters, such as the ratio of the two metric potentials. 
Currently tests of gravity on astrophysical scales are in the early stages — we 
summarize these tests and discuss the interesting prospects for new tests in the 
coming decade. 
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1. Introduction 

Einstein's theory of general relativity (GR) has proven spectacularly success- 
ful over 90 years of experimental tests [T|- These tests range from millimeter 
scale tests in the laboratory to solar system tests and consistency with gravity 
wave emission by binary pulsars. Recently cosmological motivations for mod- 
ifying gravity have renewed interest in testing GR on scales of galaxies and 
large-scale structure. 

Over the last decade it has become clear that the energy contents of the 
universe pose a major puzzle, in that GR plus the Standard Model of parti- 
cle physics can only account for about 4% of the energy density inferred from 
observations. By introducing dark matter and dark energy which account for 
the remaining 96% of the total energy budget of the universe, cosmologists have 
been able to account for a wide range of observations, from the overall expansion 
of the universe to the large scale structure of the early and late universe 2J. 
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The A-Cold Dark Matter (ACDM) paradigm has thus emerged as the stan- 
dard model of cosmology. While its phenomenological success is undeniable, 
few would argue for its aesthetic appeal. The model requires introducing exotic 
components of matter and energy. And while the vacuum energy of quantum 
fields offers a natural candidate for dark energy, the predicted amplitude from 
field theory calculations is many orders of magnitude larger than inferred from 
observations [3]. 

A critical assumption in the standard picture is the validity of GR at galactic 
and cosmological scales. Since GR has not been tested independently on these 
scales, a natural alternative is that there are new gravitational degrees of free- 
dom associated with the dark energy scale, and that these degrees of freedom 
are responsible for the observed late-time cosmic acceleration. This possibility, 
that modifications in GR at cosmological scales can replace dark matter and/or 
dark energy, has become an area of active research in recent years. See [4] and 
references therein. The driving motivation for many practitioners is of course 
the vexing problem of the cosmological constant, since any solution to this long- 
standing puzzle promises to have profound implications for the nature of cosmic 
acceleration. 

This review article discusses the prospects of testing GR on the largest scales 
through cosmological observations. Doing so in a model-independent way is un- 
fortunately not without ambiguities. Unlike solar system tests, where the mat- 
ter stress-energy tensor is known (i.e., the Sun), cosmological tests inherently 
rely on assumptions about the stress-energy composition of the universe — any 
modification to the gravitational side of Einstein's equations can equivalently 
be interpreted as a dark energy contribution to the stress-energy side [H |6l [7] . 
Throughout this paper, we assume that the dark energy component is nearly 
spatially-homogeneous. Thus strictly speaking the cosmological tests discussed 
below are in fact tests of the assumption of GR plus smooth dark energy. For- 
tunately, however, nearly all of the specific modifications discussed below also 
lead to testable deviations from GR in the solar system. 

In the rest of this introductory section, we briefly review existing constraints 
on deviations from GR from solar system and laboratory tests of gravity. In 
Seed we review some of the leading proposals for modified gravity (MG) the- 
ories, designed as an alternative to dark energy to explain the present day 
acceleration of the universe. In these models gravity at late cosmic times and 
on large scales departs from the predictions of GR. Note that for the purpose 
of this review we do not discuss alternatives to dark matter, such as Modified 
Newtonian Dynamics (MOND) |8j and its Tensor- Vector-Scalar (TeVeS) covari- 
ant generalization Olin]. In Sec. |3]we consider the prospects of testing such 
models in view of the massive new astronomical surveys and other experiments 
expected in the coming decade. We conclude with brief remarks in Sec. |4j 

1.1. Approaches to Modify Gravity 

The field equations of GR can be derived from the action 

SgR = / <i^Xy/^ R + S'mattor[.9A'i'] ' (1) 
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where Mpj = I/SttGn with Gn being Newton's gravitational constant, g is 
the determinant of the metric tensor g^^, and R is the Ricci scalar. The first 
term is the Einstein-Hilbert action, while 5'inatter contains all matter fields, with 
minimal couplings to g^^,. The stationary point of the variation of S'gr with 
respect to the metric then yields the Einstein field equations. 

The simplest and best-studied modifications to General Relativity are scalar- 
tensor theories, which introduce a scalar cousin to the graviton. First pro- 
posed by Brans and Dicke [11] as an attempt to reconcile gravity with Machian 
ideas [12l [131 IHl US] , scalar-tensor lagrangians are now understood to arise as 
low-energy limits of various theories of particle physics. These include dimen- 
sional reductions of higher-dimensional gravity and brane- world models |16l 1171 
[TSl Uni [H] . Moreover, a host of more intricate modifications of gravity, in- 
cluding massive gravity and the Dvali-Gabadadze-Porrati (DGP) model ^22l[23l 
[24l [25] , reduce to a scalar-tensor form in certain decoupling limits ^264 |27j |28j . 

The Brans-Dicke action [TI] is characterized by a universally-coupled scalar 
field. 

The effective Newton's constant is thereby promoted to a scalar field. The 
field redefinitions = ^g^u and wbd log $ map ^ to the 

Einstein frame, where the coefficient of the Einstein-Hilbert term is constant: 

^Ei^nstein ^ f ^.^^^ f^j^ _ 1 (g^)2\ ^ g^^^^^^ [g,^^^^l^'^\ . (3) 



For future use, we have introduced 2/3 = l/-\/3/2 + wbd- Unfortunately, as 
reviewed below, the Brans-Dicke parameter is so tightly constrained by solar 
system and pulsar observations, wbd ^ 4 x 10"* [T] , that the cosmological effects 
of the scalar field are rendered uninterestingly small. Since this stringent con- 
straint only applies to visible matter, many authors have relaxed the assumption 
of universal coupling and explored models in which the scalar field only sub- 
stantially interacts with the dark matter [ ^ [5g i |5T]|5 ^l551IMl [55 1 [55 1 [57 ] . 

The alternative gravity theories discussed in Sec.[2]all reduce to scalar-tensor 
theories in certain limits. (See |38j . however, for an IR modification of gravity 
without new degrees of freedom.) To produce interesting effects cosmologically 
while satisfying the tight local constraints on deviations from GR, these scalar 
fields must somehow be hidden or screened in the local environment. There are 
three known mechanisms for screening a scalar field, all of which exploit the fact 
that the matter density at a typical location in the solar system or in pulsar 
environments is many orders of magnitude larger than the mean cosmic density. 
The success of these models relies on a nonlinear mechanism that operates on 
small scales or inside dense objects to recover GR on solar-system/pulsar scales 
where gravity is well tested. The large local background density triggers non- 
linearities in the scalar field, which in turn result in its decoupling from matter. 
In short, the three mechanisms rely on giving the scalar field a large mass, a 
large inertia, or by weakening its coupling to matter, respectively. 
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The first mechanism, discussed in Sec. 2.1 arises in chameleon field theo- 
ries [3nill01llTl|32]- By adding a suitable self-interaction potential V{(j)) to ([s]), 
the chameleon scalar field acquires mass which depends on the density. The 
mass is large in regions of high density, thereby suppressing any long-range 
interactions. (Density-dependent effective couplings were initially noted in a 
different context [13] ■) Chameleon fields can in fact have different couplings 
to different matter species, thereby generalizing the universal coupling assumed 
in Brans-Dicke theories. A subclass of chameleon theories are so-called f{R) 
theories of gravity |44l|45l|46l|47l|48l|49l[50] — for ubb = and suitable choice 
of V{4>), the scalar-tensor action can be mapped through field redefinitions into 
an action of the form ([I]) with the Einstein-Hilbert term replaced by /(i?), a 
general function of the Ricci scalar. 

The second mechanism for hiding a scalar is achieved with symmetron fields 



(Sec. 2.2), proposed recently in |5ll[52]. The symmetron Lagrangian is quali- 



tatively similar to that of chameleon models, but the mechanism and its phe- 
nomenological consequences are drastically different. The screening in this case 
relies on a scalar field acquiring a vacuum expectation value (VEV) that is small 
in high-density regions and large in low-density regions. An essential ingredi- 
ent is that the coupling to matter is proportional to this VEV, so that the 
scalar couples with gravitational strength in low-density environments, but is 
decoupled and screened in regions of high density. This is achieved through a 
symmetry-breaking potential, hence the name symmetron. 

The third mechanism, discussed in Sec. 2.3 relies on the scalar field having 
derivative interactions that become large in regions of high density or in the 
vicinity of massive objects. Perturbations of the scalar in such regions acquire a 
large kinetic term and therefore decouple from matter. Thus the scalar screens 
itself and becomes invisible to experiments. This so-called Vainshtein effect [S31 
[Ml [55] ensures the phenomenological viability of brane-world modifications of 
gravity |22l(23llM1,(5^[57llMl[Ml[Mll6lllMllMllM[6g galileon scalar field 
theories [66, 67, 68, 69l|70l|711 US]- 

These various theories and the corresponding observational constraints will 
be discussed respectively in Sees. [2]and|3| To set the stage for these discussions, 
we briefiy review the current constraints on GR from laboratory and solar system 
tests of gravity. For a more detailed review of the subject, see [T]. 



1.2. Newtonian Tests 

Einstein's gravity theory is based on the weak equivalence principle (WEP), 
which states that the trajectory of a freely falling test body is independent of 
its internal structure and composition. The WEP is of course not unique to GR 
— any theory whose matter fields couple to a unique metric tensor, e.g. the 
Brans-Dicke theory in ^ , satisfies the WEP, independent of the field equations 
governing this metric. (As we review below, scalar-tensor theories, including 
Brans-Dicke, violate a stronger version of the equivalence principle satisfied in 
GR — the Strong Equivalence Principle (SEP) — which extends the universality 
of free-fall to include gravitational self-energy contributions.) 
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Eotvos-type experiments test the WEP by measuring the ij parameter — 
the fractional difference in the acceleration of freely falling bodies of different 
composition. The best limit comes from the torsion balance experiment of the 
University of Washington (Eot-Wash), which gives 77 = (0.3 ± 1.8) x 10~" [75] . 
The forthcoming space-based mission MicroSCOPE [74] is designed to reach 
10~^^ sensitivity. Other satellite experiments, such as Galileo Galilei [75] and 
the Satellite test of the Equivalence Principle (STEP) [76 , will further improve 
on current limits. 

Even when the WEP is satisfied, modifications of gravity are also constrained 
by tests of the gravitational inverse-square law [77|. The simplest example is a 
scalar field with mass m mediating a force of range A and coupling strength a, 
corresponding to the Yukawa potential 

ih = -a e ' (4) 

r 

Experimental tests can be characterized as providing limits on A and a [73l [78l 
[79I. For gravitational-strength coupling (a ^ 0{1)), there is no evidence of a 
fifth force down to A = 56 /im ,78) . 

1.3. Post- Newtonian Tests 

The early experimental successes of GR were the consistency with the peri- 
helion shift of Mercury and the deflection of light by the Sun. Beginning around 
1960, a number of experimental tests of GR were carried out, including the mea- 
surement of gravitational redshift. The consistency of the decrease in the orbital 
period of the Hulse- Taylor binary pulsar with the GR prediction of energy loss 
due to gravity waves further established the validity of GR in the weak field 
regime. 

Solar system tests of alternate theories of gravity are commonly phrased in 
the language of the Parameterized Post-Newtonian (PPN) formalism. The PPN 
formalism applies to metric theories of gravity in the weak field, slow motion 
limit, i.e. when potentials and velocities are small: ^',w^/c^ ^ ^ 1. The 
metric can then be written as a perturbation about the Minkowski metric (or, as 
we shall see, the FRW metric for the expanding universe). Terms up to second 
order in the potentials are included, with coefficients that are unity or zero for 
GR replaced by parameters that accommodate MG theories. 

The 7 and /3 PPN parameters are given by the metric 

ds^ = -(1 + 2* - 2/3*2) ^^2 + (1 _ 27*) df2 (5) 

where the potential ^E* = —G-^M/r for the Schwarzschild metric. The parameter 
7 describes the spacetime curvature induced by a unit mass and /3 the nonlin- 
earity in the superposition law of gravity. We will employ a similar metric in 
Sec. 3 with the FRW metric and an arbitrary (small) potential. In general for a 
fluid description of matter, allowing for generic Poisson-like potentials, the PPN 
formalism requires ten parameters for a complete description [1] . Note that the 
PPN formalism has constant parameters and therefore does not accommodate 
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Yukawa-like modifications with finite A. For astrophysical tests we will work 
with effective parameters that may have a scale and time dependence, but a 
similar expansion of the metric is still useful. 

The Brans-Dicke theory given by ^ has identical PPN parameter values 
as in GR, except for 7 = (1 -I- (u;bd)/(2 -I- i^bd)- Since 7 is the most relevant 
PPN parameter for the modified gravity theories of interest, we focus on this 
parameter for the rest of our discussion. 

The tightest constrain on 7 comes from time-delay measurements in the solar 
system, specifically the Doppler tracking of the Cassini spacecraft, which gives 
7 — 1 = (2.1 ± 2.3) X 10~^ [5D]. Light deflection measurements, meanwhile, 
constrain 7 at the 10"'' level [ST]. The observed perihelion shift of Mercury's 
orbit sets a weaker limit of 10"'^ 

As mentioned earlier, GR satisfies an extended version of the equivalence 
principle, the SEP, which states, in particular, that macroscopic objects follow 
the same trajectory in a uniform gravitational field as test masses. In other 
words, the universality of free-fall is preserved in GR even when accounting for 
self-gravity contributions. The SEP is violated in Brans-Dicke theories and in 
all of the MG theories considered here. Violations of the SEP result in the 
Nordtvedt effect [83] — a difference in the free-fall acceleration of the Earth 
and the Moon towards the Sun, which is detectable in Lunar Laser Ranging 
(LLR). (Because the Earth and the Moon have different compositions, however, 
one must worry about a fluke cancellation between WEP and SEP violations 
in this measurement. To disentangle these effects, laboratory tests of the WEP 
have been performed using tests masses with Earth-like and Moon-like com- 
positions [84].) Searches for the Nordtvedt effect in LLR data constrain PPN 
deviations from GR at the 10~^ level [T]. 



2. Modified Gravity Theories 

In this Section we review various proposals for modifying GR at large dis- 
tances. The list is by no means exhaustive — our goal is to provide the 
reader with an overview of a few broad classes of theories that have attracted 
considerable interest over the last few years, highlighting key theoretical and 
observational differences among them. For this purpose, we find it useful to 
group theories based on the three qualitatively different non-linear mechanisms 
through which GR is approximately recovered locally. These are the chameleon 



(Sec. 2.1), Vainshtein (Sec. 2.3) and symmetron (Sec. 2.2 1 mechanisms 



2.1. Chameleon/ f {R) Field Theories 

Chameleon field theories generalize ([3| to include a scalar potential V(0), 
whose properties will be discussed shortly, as well as a more general coupling 
A{(t)) to matter fields: 



^cha 
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Figure 1: The chameleon effective potential Vcg (solid curve) is the sum of two contributions: 
the actual potential V{(f>) (dashed curve), plus a density-dependent term from its coupling to 
matter (dotted curve). 



One can allow different couplings to the various matter fields, thereby intro- 
ducing violations of the Equivalence Principle. Furthermore, one can also in- 
troduce a coupling to the electromagnetic field strength, which induces photon- 
chameleon mixing in the presence of magnetic fields [85l [86] . For the purpose 
of this review article, we focus on the simplest case of a universal, conformal 
couphng of the chameleon, as in ([6|. 

The equation of motion for (/) that derives from this action is 

□ (/. = F^-A3(0)4^f , (7) 

where T = g^iyT^'^ is the trace of the energy-momentum tensor defined with 
respect to the Jordan- frame metric g^^, — A'^{<p)gfj_i,. Since matter fields couple 
minimally to ^^i/, this stress tensor is covariantly conserved: Wf^T^'^ = 0. 

To study the field profile on solar system and galactic scales, we can ap- 
proximate the metric in ([t]) as flat space, ignore time derivatives, and focus on 
the case of a non-relativistic, pressureless source. In terms of an energy density 
p = A^{(f>)p conserved in Einstein frame, we obtain 

V^(b = V^ + A^^p. (8) 

Thus, because of its coupling to matter fields, the scalar field is affected by 
ambient matter density. Its profile is governed by an effective potential 

Kff(0) = ^(0) + ^(</')p. (9) 
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For suitably chosen V{(j)) and A{(j)), this effective potential can develop a 
minimum at some finite field value ^min in the presence of background mat- 
ter density, as illustrated in Fig. [l] where the mass of the chameleon field is 
sufficiently large to evade local constraints: 

V<P<pi4'min) + ^,00(</'min)p ■ (10) 



TO 



2 



mm 



Assuming A{(f>) is monotonically increasing, for concreteness, the general condi- 
tions that V{(j)) must satisfy are [39l[87|: (i) to balance the potential against the 
density term, we must have V^^ < over the relevant field range; (ii) since V^^^ 
typically gives the dominant contribution to the mass term, stability requires 
V<p<p > 0; (iii) the effective mass increases with density provided that V^,^,^, < 0. 

As seen from Fig. [l| for suitable choices of V and A, the effective mass can be 
a steeply growing function of the ambient density. The chameleon mechanism 
therefore exploits the large difference between local (~ 1 g/cm'^) and cosmic 
(w 10^^° g/cm^) densities to generate a wide range in mass scales. As we 
have seen in Sec. |1.2[ tests of the inverse-square law constrain the range in the 
laboratory to be < 50 /Ltm. Meanwhile, we would ideally like to achieve rries ~ 
Ho at cosmological densities, such that the chameleon behaves as a quintessence 
field. In all known examples, however, the best one can do is rrics ^ 10'* Hq, 
corresponding to a range < 0.1 Mpc in the cosmos. The reason for this is simple 
— for a potential for which rrictf ^ Hq at cosmological densities, the field value 
where this is reached is generally a distance ^ Mp\ away from the local value, 
and maintaining gradient of such a large field difference is just energetically too 
costly. The largest gradient energy that can be tolerated corresponds in practice 
to TOcff > 10^ Hq at cosmological densities. 

A prototypical potential that satisfies all of the above conditions is the in- 
verse power-law form, 

= ' (11) 

where n is some positive constant. This falls within the class of tracker potentials 
relevant for quintessence models of dark energy [BS, ■ For the coupling function, 
a nice choice that makes contact with Brans-Dicke theories is the exponential 
form of ^: 

A{^) = e'3'^/*^P' . (12) 

The parameter /3 is implicitly assumed to be 0(1), corresponding to gravita- 
tional strength coupling. Remarkably, it was pointed out in |42| that much 
larger couplings are allowed by current constraints, but one must be concerned 
with adiabatic instabilities for /3 ^ 1 [89^. Assuming (3<j> <^ Mpi, which will be 
the case for most situations of interest, the effective minimum lies at 



PI 



/3p 

Meanwhile, the mass of small fluctuations around this value is 

n + 2 



(13) 



4+n 



TO^i„ « nin + l)M-iT^ ( j . (14) 
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Both (pmin and mmin manifestly depend on the matter density — m^i^ grows 
larger for increasing p, as desired, while (fimin decreases. 

Potentials with positive powers of the field, V^cf)) ~ 0^" with a an integer 
> 2, are also good candidates for chameleon theories. In the simplest case of 
V{(j)) = A(/>'*, [10] showed that existing laboratory constraints at the time were 
satisfied for (3 = 1 and A — 1. Moreover, the authors of [ID] identified various 
signatures for future laboratory tests of the inverse-square-law. Subsequent 
analysis by the Eot-Wash group [79] excluded a significant part of the parameter 
space. 

2.1.1. Thin-Shell Screening 

The density-dependent mass immediately results in a further decoupling 
effect outside sufficiently massive objects, due to the thin-shell effect. Consider 
a spherical source of radius R and density pin embedded in a homogeneous 
medium of density pout • The corresponding effective minima will be respectively 
denoted by (/)min-in and ^min-out- For a sufficiently massive source, the scalar 
field is oblivious to the exterior matter and is therefore pinned near 0min-in 
in the core of the object. Of course, (j) must deviate substantially from 0min-in 
near the surface of the object since (f) must eventually reach the asymptotic value 
</'inin-out far away. Thus the gradient in builds up only within a thin-shell of 
thickness Ai? below the surface, given by 

A-/? 1 f/'min— out *^min— in ^\ 

~R ^ 6/JAfpi ' ^ ' 

where $n is the surface Newtonian potential. In other words, the shell thick- 
ness is determined by the difference in (j) values relative to the difference in 
gravitational potential between the surface and infinity. 

Since field gradients are essentially confined to the shell, the exterior profile 
is suppressed by a thin-shell factor: 

\- 0min-out ■ (loj 



47rMpi R r 

The suppression factor of AR/R can alternatively be understood in a more 
intuitive way as follows. Deep inside the source, the contribution to the exterior 
profile from infinitesimal volume elements are Yukawa-suppressed due to the 
large effective chameleon mass in the core. Only the fiux lines from within the 
thin-shell can propagate nearly unsuppressed to an exterior probe. 

Clearly the thin-shell screening breaks down for sufficiently small objects. 
Imagine shrinking the source radius keeping the density fixed. Eventually, the 
cost in gradient energy to maintain the field difference between ^min-in and 
</'inin-out bccomcs too large, and the scalar field no longer reaches (?!>min-in m 
the core of the object. In this limit the thin-shell screening goes away, and the 
exterior profile takes its usual form 

' f?*min— out ■ 

(17) 



^no screen ~ * n r 

47rMpi 
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The criterion for thin-shell screening to be effective is for the right-hand side 



of (15) to be ^ 1. Clearly, the more massive the source, the easier is this 
condition satisfied, as expected. But note that the criterion also depends on 
the density contrast — for fixed source, a denser environment implies smaller 
<^min-out: which makes the thin-shell condition easier to satisfy. We will see 
in Sec. |2.1.3| that the thin-shell screening effect leads to striking observational 
signatures. 

2.1.2. Chameleon Constraints from Tests of Gravity 

The tightest constraint on the model comes from fifth-force searches in the 
laboratory [77|, which set a limit of « 50 /im [75] on the fifth-force range for 
gravitational strength coupling. Modeling the chameleon profile in the Eot- 
Wash set-up, taking into account that torsion-balance measurements are per- 
formed in vacuum, it was shown in |39j that, for the inverse power-law potential 
— M^^"^ /(f)^ with n and /? of order unity, the fifth-force constraint trans- 
lates into a bound on M : 

M < 10"^ eV . (18) 

Remarkably, the upper bound coincides with the energy scale of dark energy: 
A ~ (10-3 eV)4. This bound not only ensures that the chameleon is consistent 
with fifth force searches, but it is also sufficient to satisfy all known tests of 



GR (see Sees. |1.2| and 1.3), such as searches for WEP violation. Lunar Laser 



Ranging observations and other tests of post- Newtonian gravity [321 EO] • 

2.1.3. Observational Signatures 

Because objects of different mass have different effective coupling to the 
scalar, macroscopic violations of the Equivalence Principle can arise, even in 
the case where the microphysical lagrangian ^ displays no such violations [39l 
HOI [50] . In particular, test masses that are screened in the laboratory may 
not be screened in space. This leads to the striking prediction that satellite 
tests of gravity, such as the planned MicroSCOPE [74] , Galileo Galilei [7S] and 
STEP [TH] missions, might observe violations of the WEP with rj ^ 10"^'^, in 



blatant conflict with laboratory constraints. Meanwhile, from (17) the total 
force (gravitational -I- chameleon-mediated) between unscreened particles is a 
factor of 1 -I- 2/3^ larger than gravity, which would appear as 0(1) deviations 
from Gn measured on Earth. 

Effective violations of the WEP also results in a host of astrophysical sig- 
natures uncovered in [90 . For example, consider (unscreened) dwarf galaxies 
in large void regions. The HI gas in these galaxies is also unscreened, but the 
stars are screened. This should result in a systematic 0(1) mismatch in the 
rotational velocities of these two tracers, and hence a corresponding mismatch 
in mass estimates. 

Other interesting observable signals arise if the chameleon couples to the 
electromagnetic field 

J d'^x^e^-'I'/^'^^'F^^F'"' . (19) 
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Such a coupling results in photon-chameleon oscillations in the presence of an 
external magnetic field. The GammeV experiment (HU E2 ES] searches for the 
afterglow [Ml [23 that would result from trapped chameleons converting back 
into photons and has so far excluded the range 5 x 10^^ < P < 6.4 x 10^^ W^ - 
This bound has been further refined recently in a more complete analysis [55] . 
The second-generation experiment, GammeV-CHASE, is currently taking data 
and will further improve on this bound |93) . Photon-chameleon mixing can 
also be detected astrophysically, through induced polarization in the spectrum 
of astrophysical objects [96] and enhanced scatter in the X/7-ray luminosity 
relation of active galactic nuclei (AGN) [37]. A recent analysis of AGN data 
shows tantalizing evidence of this effect [97j . 

2.1. 4- Cosmology of Chameleon Models 

The chameleon scalar field can act as dark energy driving cosmic accelera- 



tion. Unfortunately the bound ( 18 ) from local tests of gravity forces the equation 
of state of the chameleon component to be indistinguishable from w = —1, re- 
sulting in an expansion history identical to ACDM. This is closely related to a 
point made earlier: in order to generate observable deviations from w — ~1, the 
chameleon field must be rolling significantly on a Hubble time, which requires 
TOeff ~ Hq on cosmological scales. But this in turn requires a field displacement 
of order Mpi from the local value to the cosmos, which is too costly energetically 
— galaxies can no longer maintain their thin shell, leading to gross violations 
of local constraints. 

Explicitly, the inverse power-law form, V{(p) — A/^+"/(/i", leads to too small 
a dark energy contribution for M ~ 10"'^ eV. In order to account for cosmic 
acceleration, we must lift this potential by an appropriate constant, which of 
course has no effect on our earlier analysis of local tests. We can exploit the 



coincidence of the dark energy scale and the upper bound in ( 18 1 to choose a 
potential with a single mass parameter |41j 

= M4exp(M"/(/)"), (20) 

with M = 10-3 eV. For > M, this gives V{(j)) « M"* -I- M^+'' -f . . ., 



which is (11) plus a constant. Thus, while cosmic acceleration is not explained 
in chameleon models, it is easily incorporated within its phenomenology. 

Let us now turn to the cosmological evolution of the chameleon. In the 
matter-dominated era, ([t]) reduces to 

4> + 3H<j) = - A,4,p„, , (21) 

where satisfies the usual conservation law: pm + SHp^^ ~ 0. In particular, 
the effective potential for is time-dependent. The behavior of the scalar can 
be inferred by standard adiabaticity arguments — as long as mmin 3> H ^ the 
field continuously adjusts to the evolving minimum. This is manifestly satisfied 



for the potential (20 1 and coupling function (12 1, 



~3/3r!„,^<il + n + n( » 1 , (22) 
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since cfiniin ^ -^^Pi for relevant densities, as mentioned earlier. The fact that 
w((/)niin) 3> -ff on cosmological scales today is yet another way of seeing that 
chameleon dark energy is indistinguishable from a cosmological constant at the 
background level. 

The story is more subtle in the radiation-dominated era, since w for 
a relativistic fluid, and hence the source term in ([t]) vanishes. This suggests 
that the field is overdamped and remains frozen until matter-radiation equal- 
ity, which would lead to inconsistencies with nucleosynthesis if the field started 
out in the early universe with generic initial displacement ~ Mp\. Fortunately, 
the trace is not always small for a realistic relativistic fluid. As the universe 
expands and cools, matter species in the thermal bath successively become non- 
relativistic when m ^ T. Whenever this happens, becomes non-zero for 
about one e-fold of expansion, until the species in question becomes Boltzmann 
suppressed, thus driving the field closer to the minimum [98]. The total dis- 
placement incurred by "kicks" from known particles until nucleosynthesis is of 
order (3M-p\ |3T]. Thus, for a broad range of initial conditions the field settles 
to the effective minimum by the onset of nucleosynthesis. 

2.1.5. f{R) Theories 

A class of chameleon theories that has sparked considerable interest over the 
last few years is the f{R) form [HJHS], where the Einstein-Hilbert term in ([I]) 
is replaced by a general function of the Ricci scalar [44l |45l |46] : 

^#2 f 

5 = ^ y d^x^f{R) + 5„,attcr[5M.] • (23) 

See [47l |48] for recent reviews. This action was first introduced by Starobinsky 
in the context of inflation [99] and studied further in [1001 IIOH I102j . Functions 
of other curvature invariants have also been considered [103] , but this generally 
introduces ghost instabilities [104) . (Of course general higher-order curvature 
invariants are harmless from an effective field theory point of view provided they 
contribute small corrections to Einstein gravity. But here one is precisely inter- 
ested in order-unity deviations from General Relativity.) The only exceptions 



are functions of the Ricci scalar (as in (23)) or functions of the Gauss-Bonnet 
term [TDS]. 



It is well-known that (23 1 is just a scalar-tensor theory in another guise. To 



see this, introduce a non-dynamical scalar ip: 

S^^j d'^V^ [.fW + %(^-^)]+ 5mattcr[5^.] ■ (24) 

Varying with respect to ip imposes the constraint ip — R, assuming /" 7^ 0. And 
since the constraint is algebraic, we can substitute it back into ([24]), thereby 



reproducing (23). Furthermore, the field redefinitions Qf^i, — {df /d'ip)gfj,i, and 



= — ■\/3/2Mpi logd//dV' reduce (24 1 to the more familiar form 



d'xV^ { ^i? - ^{dcj^r - y(0) ) + ^„.attcr g^.eV2/30/Mp, ^ (25) 
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where the only remnant of the f{R) choice is encoded in the scalar potential 



V 



\d4, J 



(26) 



The action ([25| is precisely of the form Theories of f{R) gravity are there- 
fore physically equival ent to a subclass of chameleon field theories, namely those 
with A{(j)) given by (12) for the specific choice (3 — 1/a/6, as can be seen 
from (25). In other words, only for this specific choice of /3 do chameleon 



theories allow an f{R) description. 

For completeness, we give the equations of motion in the f{R) language. 
Varying ( 24 ) with respect to the metric gives the modified Einstein field equa- 
tions 

Gp. + f,RR,.u - ({ - .9^- - ^^.^uf,R = SnG^Tf,, , (27) 

where T^^ = —{2/^J—g)5S-aiattcil5g^^ is the matter stress-energy tensor, and 



= df/dR. Taking the trace of (27) yields an equation for the scalar degree 
of freedom f^n: 



(28) 



Two popular choices of f{R) are the Hu-Sawicki |49| and Starobinsky [50] 
forms, given respectively by 



f{R) 
f(.R) 



R 



(HS) 



R + aAV 



1 



-a/2 



- 1 



(29) 
(30) 



where a and a are positive constants. (For convenience we have redefined the 
parameters of |3S] as follows: M"^ = ni^/c. 



l/r 



-1 + 1/n 
C1C2 



Both forms 



satisfy /(i?) — )■ as i? — )■ 0. For the relevant choices of parameters, however, 
the entire cosmological history of our universe lies in the R 3> Af ^ regime, in 
which both forms have identical limits: 



f{R)KR + aM^ 



R 

M2 



{R > Af2) . 



(31) 



In particular, this is the relevant regime to describe the chameleon mechanism. 
Using the above dictionary, the corresponding potential for the canonically- 
normalized scalar field in Einstein frame is 



F(0) 




(32) 
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This potential is finite at </) = (corresponding to i? oo in the f{R) language), 
but its derivative blows up at this point for any a > 0. This is a generic 
consequence of starting with a nice functional form in the f{R) language — 
translating to the scalar-tensor language involves inverting a relation to obtain 
R = R{4>) , which therefore typically results in a non-analytic scalar potential. 

The pathological form of the potential is indicative that quantum correc- 
tions become important near = [MQ The scalar potential V{<l)) receives 
radiative corrections from (j) loops, whose finite part at 1-loop is 

A..-.,^S|.o.(^)^(^)-=«^,„.(^), ,33, 

(As usual there are also divergent pieces that must be fine-tuned to zero.) This 
overwhelms the classical part for (p <C Mpi. Thus scalar configurations in the 
small-field limit are in the quantum regime in this case. 

The cuspy form of the potential leads to pathologies already at the classical 
level |1071 1108[ I109j : although the static chameleon solutions described earlier 
only rely on becoming large, the fact that V itself remains finite aX (j) — 
implies that finite-energy configurations {e.g. a collapsing star) can dynamically 
overshoot the effective minimum and hit the singular point. The singularity at 
= also led [llOj to argue against the existence of relativistic stars in f{R) 
gravity, though this conclusion has been questioned recently by . What p3| 
shows, however, is that the classical approximation is not valid for <^ Mpi — 
the f{R) dynamics are dominated by quantum effects in this regime. 

We should stress that these problems are specific to cuspy potentials of 



the form ( |32[ ). The inverse power-law form of (111, for instance, is clearly 
immune to singularity issues. If a dynamical configuration causes the field to 
overshoot, the infinite barrier at = will stop the field and lead to oscillations 
around the minimum. (Alternatively, in |112[ 1113] it was shown that adding 



higher-order curvature invariants to (31) can remove the singular point and 



push the high-curvature regime to infinity in field space.) Similarly, the 1- 
loop contribution ( |33[ ) can be under control over the relevant field range in this 
case |106j . 

What choice of f{R) leads to inverse power- law potential? In the strong 
curvature regime, the desired form is 



f{R) = R- aM^ 



{R > M^) . (34) 



This is similar to (31), except that the power is positive in this case. Moreover, 
whereas (31) can be an analytic function of R, the above has branch cuts for 
any choice of n. As mentioned earlier, this is an unavoidable consequence of 
translating between the f{R) and scalar-tensor forms — analytic V^(0)'s map 



^We are grateful to Andrew Tolley for helpful discussions on this point. 



15 



onto non-analytic /(i?)'s, and vice versa. The stability discussion above there- 
fore implies that non-analytic dependence on the field variables is forced upon 
us if we insist on working in the f(R) description. 

2.2. Symmetron Fields 

A second mechanism for hiding scalar fields was proposed recently in the 
context of symmetron field theories |511 152) . As we will shortly, the symmetron 
technology and building blocks are similar to those of chameleon models. But 
the physics of the screening mechanism and its phenomenological consequences 
are dramatically different. In particular, unlike in chameleon theories, the sym- 
metron has a small mass everywhere and therefore mediates a long-range force 
on Earth and in the solar system. 

In the symmetron mechanism, the vacuum expectation value (VEV) of the 
scalar depends on the local mass density, becoming large in regions of low mass 
density, and small in regions of high mass density. In addition, the coupling of 
the scalar to matter is proportional to the VEV, so that the scalar couples with 
gravitational strength in regions of low density, but is decoupled and screened 
in regions of high density. 

The starting point is the general chameleon action ([6]). In particular, we 
assume that the symmetron couples universally to matter fields through a con- 
formal coupling. The scalar potential is of the symmetry-breaking form, 

V{cf>) = ~^fi'4>^ + ^Xcf>\ (35) 

involving a tachyonic mass scale fi and dimensionless coupling A. Meanwhile, 
the coupling to matter is governed by 

Aicb) = l+^^'+0{<j>'/M^). (36) 
The field range of interest is restricted to (/> ^ M, as we will see shortly, thus 



higher-order terms in A{(j>) are negligible. The theory described by ( 35 ) and ( 36 1 
has a 1i2 symmetry -> —(p. The self-interaction potential is the most general 
renormalizable form consistent with this symmetry. The quadratic coupling to 
the matter stress tensor is the leading such coupling compatible with the I12 
symmetry. It is non-renormalizable, suppressed by the mass scale M. 

The screening mechanism is achieved through an interplay between the 
symmetry-breaking potential and the universal coupling to matter. In the case 
of a non-relativistic source, the effective potential ^ is given by 

VM<^) = l{J^~^,')c^' + \\c^\ (37) 

Thus, whether the quadratic term is negative or not, and hence whether the Z2 
symmetry is spontaneously broken or not, depends on the local matter density. 
In vacuum, the potential breaks the reflection symmetry spontaneously, and the 
scalar acquires a vacuum expectation value (VEV) 4>o = fi/VX. Inside a source 
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of sufficient density, such that p > M^fj,'^, the effective potential no longer breaks 
the symmetry, and the VEV goes to zero. 

An essential feature is that the lowest order coupling of matter to the sym- 
metron is ~ pcfi^ /M"^ . Fluctuations 5(j) around a local background value 0veV: 
as would be detected by local experiments, therefore couple as 

^^^^'^P' (38) 

that is, the coupling is proportional to the local VEV. In high-density, symmetry- 
restoring environments, the VEV should be near zero, and fluctuations of (j) 
should not couple to matter. In rarified environments, where p < AI^p^, the 
symmetry is broken and the coupling turns back on. 

The case of interest is where the field becomes tachyonic around the current 
cosmic density: H^M^y ~ p^M"^. This fixes p in terms of M, and hence the 
mass mo of small fluctuations around 0o = p/VX- 

mo^V2p^ . (39) 

Local tests of gravity, as we will see, require M < lO^^Afpi. Hence the range 
ttiq ^ of the symmetron-mediated force in vacuum is < Mpc. Meanwhile, if this 
extra force is to be comparable in strength to the gravitational force in vacuum. 



then from (381 we must impose 0o/Af^ ~ 1/A/pi, that is, 

(40) 



p M2 

7^^ Mpi 



Together with (|39|, this imphes A - M^^H^/M^ < 1. (For M = IQ-^Afpi, 
in particular, this gives A — 10^^^^ — the symmetron is extremely weakly 



coupled.) We see immediately from (401 that 0o ^ hence the field range 
of interest lies within the regime of the effective field theory, and 0{4>'^/M'^) 



corrections in (36) can be consistently neglected. 



2.2.1. Symmetron Thin- Shell Screening 

Symmetron solutions around a source display a thin-shell effect analogous 
to the chameleon behavior discussed in Sec. |2.1.1[ Consider once again the 
ideal case of a static, spherically-symmetric source of homogeneous density p > 
p^M"^ . For simplicity, we further assume the object lies in vacuum, so that the 
symmetron tends to its symmetry-breaking VEV far away: — 0o as r — )■ oo. 

For a sufficiently massive source, in a sense that will be made precise shortly, 
the solution has the following qualitative behavior. Deep in the core of the ob- 
ject, the symmetron is weakly coupled to matter, since the matter density forces 
« there. Near the surface, meanwhile, the field must grow away from (j) — 
in order to asymptote to the symmetry-breaking VEV far away. The symmetron 
is thus weakly coupled to the core of the object, and its exterior profile is dom- 
inated by the surface contribution. In other words, analogously to chameleon 
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models, there is a thin shell screening effect suppressing the symmctron force 
on an external test particle. 

The parameter that determines whether a solution is screened or not is 

(Recall that <i>N is the surface gravitational potential.) Objects with a ^ 1 
display thin-shell screening, and the resulting symmetron-mediated force on a 
test particle is suppressed by 1/a compared to the gravitational force. Objects 
with a <C 1, on the other hand, do not have a thin shell — the symmetron gives 
an 0(1) correction to the gravitational attraction in this case. 

2.2.2. Tests of Gravity and Observational Signatures 

Since the symmetron-mediated force is long-range in all situations of interest, 
and because the symmetron couples to matter universally, the relevant tests of 
gravity are the same that constrain standard BD theories: solar system and 
binary pulsar observations. A necessary condition is for the Milky Way galaxy 



to be screened: ac ^ 1- Since $ ~ 10 for the galaxy, (41 1 implies 



M < lO^^Mpi . (42) 

It turns out that this condition is also sufficient to satisfy all current constraints. 
Indeed, as shown in j51j, with M — 10^'^Afpi the symmetron predictions for 
time-delay, light-deflection, perihelion precession of Mercury and Nordvedt ef- 
fect are comparable to current sensitivity levels and therefore detectable by 
next-generation experiments. Note that pushing M to larger values is also de- 
sirable cosmologically, since the range of the symmetron force grows with M. 



In particular, from (|39|), ^"^ < Mpc for M < 10~^Mpi. 

The symmetron is observationally distinguishable from other screening mech- 
anisms. In chameleon theory, as discussed in Sec. |2.1.2t the tightest constraint 
comes from laboratory tests of the inverse square law. Once this is satisifed, 
however, the predicted solar system deviations are unobservably small. In con- 
trast, as mentioned above the symmetron predictions for solar system tests are 
just below current constraints. On the other hand, chameleon and symmetron 
models have in common the prediction of macroscopic violations of the WEP, 
which can show up in various astrophysical observations |90j . In the Vainshtein 
case, as mentioned in Sec. |2.3.2[ the DGP model predicts modifications to the 
Moon's orbit that are within reach of next generation Lunar Laser Ranging 
observations, but light-deflection and time-delay signals are negligible. 



2.3. Theories of Massive/Resonance Gravity 

A class of infrared modifled theories that has spurred a lot of activity in- 
volves giving the graviton a small mass or width. This is alluring from a purely 
theoretical perspective, because of the formidable challenge devising consistent, 
ghost-free theories of massive gravity has proven to be. We know how to give 
a mass to all particles, but doing so consistently for the graviton remains an 
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open problem. Besides the theoretical appeal, the main motivation remains of 
course the cosmological constant problem. The cosmological constant is the 
zero-momentum component of the stress-energy tensor, hence its backreaction 
depends sensitively on the nature of gravity in the far infrared. 

At the linearized level, all theories of massive/resonance gravity reduce ef- 
fectively to scalar-tensor theories, even in the limit of vanishing mass, where the 
scalar is the longitudinal mode of the graviton. This is the famous van Dam- 
Veltman-Zakharov (vDVZ) [114] discontinuity. The resolution is the Vainshtein 
mechanism |53j — non-linearities in the longitudinal dominate in the vicinity 
of astrophysical sources, and result in its decoupling from matter. Thus GR is 
approximately recovered in the solar system. 

2.3.1. Fiertz-Pauli Gravity and Its Discontents 

The perennial challenge for theories of massive/resonance gravity is the 
avoidance of ghost instabilities. A theory of a single, Lorentz-invariant mas- 
sive graviton |115j is generically plagued with ghosts. (The story is different for 
Lorentz-violating massive gravity |1161lll7llll8| : for simplicity, we focus here on 
Lorentz-invariant mass terms.) At the linearized level, the Fiertz-Pauli tensor 



structure [115j is chosen precisely to remove the 6th polarization mode of the 
graviton, which would otherwise be a ghost. But since this is no symmetry to 
prevent this mode from being excited non-linearly, one generically finds a ghost 
propagating around non-trivial backgrounds. The instabilities of massive grav- 
ity were first diagnosed by Boulware and Deser |119| in the ADM formalism: 
integrating out the lapse function and shift vector results in a Hamiltonian that 
is unbounded from below. In recent developments, a realization of massive grav- 
ity using an auxiliary extra dimension |631 164j has been argued to have positive 
definite Hamiltonian. 

The origin of the ghost is most transparent in the Goldstone-Stiickelberg 
language |120l I12l[ I122j . As in the previous subsection, the metric fluctuation 
/i/jy = ff/ii^ ~ Vfj.i' is completed to a gauge-invariant object /i^^ by including the 
Goldstone modes. Focusing on the the longitudinal or helicity-0 mode tt, we 
have [55] 

TT 2r^ 

hfj.,y = hf,„ + -n^,^---- + -^d^,dy-K -t- YTT^t^^^'^ dydjTT , (43) 

iVjfpi Mpi -''^Pl 

where the rj^iyir/Mpi term is introduced for convenience to diagonalizc the kinetic 
matrix of h and tt. (This is the linearized version of a conformal transformation 
to Einstein frame.) Since /i^ji/ is gauge-invariant, we are allowed to add £mass = 
\/~35'"'5"^(o/i^iq/ii//3 + bh^yhap) to the Einstein-Hilbert lagrangian, where a 



and h are constants. Substituting the decomposition (43), we risk generating a 
(□tt)^ term in the action from the square of the d^d^i: term, which would yield 
a ghost at the quadratic level. But the Fiertz-Pauli structure (a = —h) is chosen 
precisely to remove this dangerous term. 

We can zoom in on the Goldstone sector by taking the decoupling limit : 
Afpi — > oo and — t- oo, keeping the strong coupling scale (r~^Mpi)^/^ fixed. 
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The resulting action PS I IT ^ . 

Sfp = J d^x (^S^Dtt + ^Dt^ [(°^)' - (d^d^n)^] + j^^^T^ , (44) 

is clearly well-defined at the quadratic level, as advocated. However, the theory 
has non-linear instabilities: expanding the interaction term to quadratic order 
in perturbations ip around a background solution n(a;) yields higher-derivative 
terms of the form □n(n(^)^, which signals a ghost. One could hope to remove 
the ghost by adding suitable higher-order terms in h can remove the ghost, 
but |120j argued there is insufficient freedom to do so. However, these conclu- 
sions have been questioned recently by |63j in the context of auxiliary extra 
dimensions. Another proposal for a non-linear version of massive gravity was 
proposed recently in |123| . 

2.3.2. Dvali-Gabadadze-Porrati Gravity 

The DGP model [22l [23l [24] shares many properties of massive gravity, but 
avoids some of its pitfalls. See f25' for a review. In this construction, our visible 
universe is confined to a 3-brane embedded in an empty, 4-Hl-diniensional space- 
time, where the extra dimension has infinite- volume. The hallmark of DGP 
gravity is the inclusion of an Einstein-Hilbert term on the brane. The action, 
given by 

-.95i?5 + / d^a^V^gl (Ra + /:mattor) , (45) 




2 

therefore involves two Planck scale, Mpi and M5. The R4 term on the brane 
gives a large inertia to the graviton and leads to a recovery of a force law 
at short distances. But since the extra dimension is infinite in extent, GR is 
not recovered in the infrared — instead the gravitational force law asymptotes 
to at large distances, corresponding to 5D gravity. The cross-over scale Tc 
from to is set by the bulk and brane Planck scales: 

Hence a large hierarchy between M5 and A/pi is necessary to recover the inverse- 
square law on sufficiently large distances. Setting Vc ^ Hq^, for instance, 
requires M5 100 MeV. 

From the point of view of a brane observer, the 5 helicity states of the 
massless 5D graviton combine to form a massive spin-2 representation in 4D. In 
the weak-field limit, the metric fiuctuation on the brane for a given an energy- 
momentum source is [22l 1124) 




h^^ = — ^V^= T^u - -v,.uT . (47) 



The infrared modification is manifest in the form of the propagator, l/(— □ + 
r""'"-^/— □), which describes a resonance graviton with tiny width r^^. (In the 
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parametrization of (69), this corresponds to a = 1/2.) As in massive grav- 
ity, the tensor structure of the above ampHtude is different than predicted by 
General Relativity, even in the Umit Tc -> oo. This is the famous vDVZ dis- 
continuity |114j of massive gravity. As we will review shortly, this discrepancy 
is an artifact of perturbation theory — non-linear effects in the helicity-0 mode 
near astrophysical sources lead to a recovery of General Relativity. This effect, 
first conjectured by Vainshtein in the context of massive gravity [S31[S11[SS], has 
been shown explicitly in DGP [Eg [1261 IMl [EH] . 

Since the brane is a codimension-one object, the standard Israel junction 
conditions lead to a Friedmann equation which is local on the brane |129) : 



H 



where the ± branches correspond to different brane embeddings in the Minkowski 
bulk space-time. The branch has generated a lot interest because it leads to 
cosmic acceleration at late times without vacuum energy [12911130] . However, it 
has been established that this branch of solutions suffers from ghost instabilities, 
first using perturbative arguments [27l [28l I131j and, more recently, by studying 
non-linear solutions |132[ 11331 11341 [T^ . The "— " branch, however, is stable. 
(To account for cosmic acceleration, the — branch requires a vacuum energy 



component. Together with the H/rc correction in (48 1, this leads to an effective 
dark energy component with w < —1 [136].) This is a generic property in DGP: 
solutions come in pairs, with one member being continuously connected to the 
trivial solution and having stable perturbations, while the other being connected 
to the self-accelerated cosmological solution and having unstable perturbations. 

As in massive gravity, we can zoom in on the helicity-0 sector via a decou- 
pling hmit [371 [SI]: -^Pi — ^ oo and Tc — )• c«, keeping the strong coupling scale 
(r-2Mpi)i/3 fixed. (See [l37] for an argument against the validity of taking a 
decoupling limit in the DGP case.) In this limit, the theory becomes local on 
the brane and describes a scalar field tt: 

Sdgp = [ d^x (snOn - -^(9^)2n7r + ^ttt) . (49) 

This theory enjoys a Galilean shift symmetry [271 [35] . 

df_,TT a^TT + , (50) 
which is a vestige of the full 5D Lorentz transformations |138j . Thus, tt has been 



dubbed a galileon field [Ml [671 [68l [Ml [iQl [711 [72] . Even though (49) contains 



higher-derivative interactions, the resulting equation of motion for tt, 

9'' (6Mpia^7r + 2r2a^7rn^ - rld^idnf) = -T , (51) 

is nevertheless second-order — all higher-derivative terms cancel out when per- 
forming the variation. This is unlike the case of massive gravity, where the 



variation of (44) yields a higher-order equation of motion. The 7r-action (51) 
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was generalized in to include higher-derivative interactions that preserve 
the Galilean shift symmetry and yield second-order equations of motion. 

The ttT coupling generates a scalar contribution to the graviton exchange 
amplitude, which is responsible for the discrepancy between the 1/3 coefficient 



in (47) and the expected factor of 1/2 for a massless spin-2 field. As mentioned 
earlier, this discrepancy persists even as Tc — t- oo. In the presence of a large 
source, however, the self-interactions terms are important and result in the 
decoupling of tt from the source. This Vainshtein screening can be immediately 



understood from (51) by considering the static, spherically-symmetric galileon 




profile due to a point source. With trivial asymptotic conditions, (51) can be 
integrated to give 

A /TV Q 1\ /f _ -/.-./ / /I (V. 3 \ 

(52) 

V tJ I f 

where we have introduced the scale, 

r. = (r^rsch)^/' , (53) 

with rsch denoting the Schwarzschild radius of the source. At short distances, 
r <^ r^, the galileon- mediated force is manifestly suppressed compared to the 
gravity: 



IWI 



< 1 . (54) 

Thus, as advocated, the strong interactions of tt lead to its decoupling near a 
source, and the theory reduces to Newtonian gravity. The deviation from stan- 
dard gravity, albeit small in the solar system, is nevertheless constrained by 
Lunar Laser Ranging observations: Tc > 120 Mpc |139( 1122] I140j . A compara- 
ble bound on has also been obtained by studying the effect on planetary or- 
bits |141j . The next generation Apache Point Observatory Lunar Laser-ranging 
Operation will improve this bound by an order of magnitude |142j . thereby 
probing the interesting regime ~ Fl^^ . 

At large distances, r 3> "r*, on the other hand, the non-linear terms in tt are 
negligible, and the resulting correction to Newtonian gravity is of order unity: 



F^ 



J- grav 



= \ ■ (55) 



The galileon-mediated force therefore leads to an enhancement of the gravita- 
tional attraction by a factor of 4/3. In this far-field regime, the theory reduces 



to a Brans-Dicke theory with wbd = 0, consistent with (47) 



Unlike chameleon theories, there are no macroscopic violations of the WEP in 



this theory [OT. This is because (51 ) takes the form of a total divergence 



-T, which follows from the shift symmetry (50). This implies a generalized 



"Gauss' law" [281 [90]: spherically-symmetric exterior solutions for tt only depend 



on the mass enclosed, as can be seen from (52). 
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The Vainshtein screening can also be understood by studying perturbations 



(/5 around a background 11 profile: n = H + ip. Expanding (49) to quadratic 
order in (f gives 

Spert = - / (3?7^, + 2K^, - 2t^^,K) d^^ipd"^ , (56) 

where K^^^, = — r^9^9j,n/Mpi is the extrinsic curvature of the brane in the de- 
coupling limit. In the non-linear regime, X S> 1, galileon perturbations therefore 
acquire a large kinetic term and decouple from matter. 

We note in passing that evaluating the above kinetic matrix for the spherically- 



symmetric background (52 1 shows that perturbations propagate superluminally 
in the radial direction in the range rgch <C r <C r^, [138] . Superluminal propaga- 
tion is a generic property of galileon theories [BS] and has recently been shown 
to hold around exact solutions of DGP |143) . The existence of superluminality 
in these models casts doubts on whether they can be completed in a quantum 
field theory or perturbative string theory |138j . 

As the above discussion illustrates, much of the DGP phenomenology (and 
its predictions for tests of gravity) is captured by the decoupling theory. This 
motivated the authors to 68J to propose a AD theory of modified gravity, by 



promoting (49) into a fully covariant, non-linear theory of gravity coupled to a 
galileon field. See also ^67 for covariantizations of more general galileon theories. 
While the extension is by no means unique, a canonical choice is 

S = J d'xV^g (^^e-2-/^^p.i? - {^^idnfDn + £„,attcr[5]) • (57) 



(The ttDtt term in (49) arises by diagonalizing the kinetic matrix as in (43), 
which also generates the ttT coupling.) Remarkably, the 4D cosmology that fol- 
lows from this action reproduces many features of the full-fledged DGP model j68] . 
The Friedmann equation has two branches of solutions, with stable or unstable 
perturbations, depending on the sign of tt. Moreover, there is a cosmological 
analogue of the Vainshtein effect: at early times, when the density of the uni- 
verse is high, non-linear interactions in tt are important, resulting in the galileon 
energy density being subdominant compared to the matter or radiation fluid. 



2.3.3. Degravitation 

One of the underlying motivations for massive/resonance gravity is the cos- 
mological constant problem. This is usually stated as "Why is the vacuum en- 
ergy so small?" But since the only observable effect of the cosmological constant 
is through its backreaction on the expansion of the universe, a more conservative 
question is "Why does the vacuum energy gravitate so weakly?" If gravity acts 
as a high-pass filter, in particular, a large cosmological constant can result in a 
small expansion rate [1441 11451 1122| . The vacuum energy can thus degravitate. 
As we will see below, degravitation is closely tied to massive gravity — any the- 
ory that exhibits degravitation must reduce, at the linearized level, to a theory 
of massive/resonance gravity |122j . Degravitation is the gravitational analogue 
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to the well-known screening of charges in a superconducting (or Higgs) phase of 
electromagnetism [122j . 

A phenomenological modification to Einstein's equations that encapsulates 
degravitation is 

G^HOrl)G^, = 8TrT^,, (58) 

where Newton's constant, G-[^{Drl), has been promoted to a high-pass filter with 
cut-off scale re- sources with characteristic wavelength <^ Vc gravitate normally, 
but those with wavelength Tc are degravitated. 



Equation ( 58 ) manifestly violates the Bianchi identity — an immediate con- 
sequence of general covariance — , and hence cannot be the whole story. But 
there is already a problem at the linearized level, g^n = rj^v -f h^^, where the 
Bianchi identity is trivially satisfied. Rewriting the filter function as Gj^^ (Dr^) = 



1 — m (□r^)/n and choosing de Bonder gauge, d^h^j^^ = di,h/2, (58) reduces to 



(□ - m^{arl)) ( V - oVf^i^^i = 8ttT^, . (59) 



Suppose for a moment that m^(nr^)) = is constant. Then the problem is 
obvious: the only allowed (Lorentz-invariant) tensor structure of a spin-2 mass 
term that is ghost-free is the Ficrtz-Pauli form |115j : mn^^h^u — hr]^^). The mass 



term in ( 59 1 is manisfestly not in this form, hence it cannot describe a consistent 
theory of a massive spin-2 particle. Allowing for more general r7i^(nr^)) only 
exacerbates things, for we can spectrally represent the graviton propagator, 

' - '~^M^?^^. (60, 



□ -m2(nr2) 7o D-M^ 

and apply the argument to each massive graviton in the continuum. 

The resolution of this paradox is simple: a massive graviton has 5 polariza- 



tion states (2 helicity-2, 2 helicity-1 and 1 helicity-0 states), but (59 1 is an effec- 
tive equation that only describes the helicity-2 part of the graviton, obtained by 
integrating out the other 3 polarization degrees of freedom. For completeness, 
let us briefly review the proof of |122( I146j . The starting point is a generaliza- 
tion of Fiertz-Pauli theory [115] , allowing for a momentum-dependent graviton 
mass: 

{^h)^^ + iK. - V^..h) = SttT^, , (61) 

where (Sh)^^ — + . . . is the linearized Einstein tensor. The diffeor- 

morphism gauge symmetry can be made explicit by introducing a Stiickelberg 
fleld [2g 

hf,u = hf,^ + df,A^ + d^Af, , (62) 

such that 6hfj,i, = df_i£_i, + d^(,fj, and SAf^ = — under gauge transformations. 
Thus hfj^i, is gauge invariant. This is analogous to electromagnetism in a Higgs 
(or superconducting) phase in which the photon becomes a gauge-invariant ob- 
servable: A^ — Afj_ + dfjcf). 
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The idea is to solve for and substitute the result back into (61 ) to obtain 
an equation for the helicity-2 modes. First, we substitute the decomposition (62 1 
into Jell): 



{£h)^, + m\Urt) [h^, - i^^.h + d^A, + d^A^ - 2r;^,c'"A„ j = SttT^, . (63) 
Taking the divergence of this expression allows us to isolate an equation for A^: 



-9'^ h 



(64) 



where = df^A^—d^A^. Note that taking another divergence yields: d^^d'^h^^- 
Dh = 0, hence h can be represented in the form h^^^ = h^i, — rjf^^Ilaph"^ /3, 
where IIq^ = ijap — dadp/O is the transverse projector. Thus h^i, carries two 
degrees of freedom. Coming back to ( 64 1 , the solution for A^, is 



A„ 



(65) 



where Q is an arbitrary gauge function. Substituting ( 65 ) back into (|63| , and 
choosing Q appropriately, we find 



1 



□ 



(66) 



which is the linearized version of (58 1. Therefore, as advocated, (59) is an 
effective equation for the helicity-2 part /i^^, after integrating out the extra 
helicities of a massive spin-2 representation. 

Hence, any theory that proposes to filter out the cosmological constant must 
reduce, in the weak-field limit, to a theory of massive or resonance gravity. Of 
course the reverse is not true — not every theory of massive/resonance gravity 
will exhibit degravitation at the non-linear level. To see degravitation at work 
at the linearized level, consider (61) for vacuum energy, T^j^ = —Arj^^, and let 
us focus on the case of constant nr for simplicity: 



(67) 



Note that without the mass term the solution grows unbounded, hij ~ A{t'^5ij + 
XiXj)/6, which is the weak-field version of de Sitter space. With the mass term 
turned on, however, the solution is just flat space 



Ar? 



(68) 



Thus the gravitational backreaction of A vanishes in a theory of a massive spin-2 
particle. 

What functional forms for m^(nr^) are allowed? A useful parametrization 
is the power-law form [1221 1146j 



2a - ^-2(i-q) 



(69) 
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where a is a constant. First, in order for the modification to be relevant in the 
infrared, we must impose a < 1, otherwise the mass term becomes neghgible as 
□ — ^ 0. Second, in order for the graviton propagator to have a positive definite 



spectral density, p(M^) > in (60 1, we must require a > 0. Otherwise the 
left-hand side of (60) vanishes as □ — )■ 0, which is inconsistent with p(M^) > 
0. Finally, in |122j it was argued that a < 1/2 is necessary in order for the 
degravitation phenomenon to be effective in a certain decoupling limit of the 
theory. Hence, the allowed range is 

< a < ^ . (70) 

The lower bound corresponds to a hard mass for the graviton, discussed in 
Sec. |2.3.1[ The range a > corresponds to a continuum of massive graviton 
states, which immediately points to extra dimensions of infinite extent. The 



best known example is the DGP model (Sec. 2.3.2), with D = 5 bulk space- 
time dimensions. The DGP mass term, m^(nr^) = r^^-^/— □, corresponds to 



a = 1/2, which coincides with the upper bound of (70). Cascading gravity 



theories (Sec. 2.3.4) generalize the DGP scenario to higher dimensions. As we 



will see in Sec. 2.3.4[ all such theories correspond to a « in the far infrared. 



2.3.4- Cascading Gravity 

It is natural to ask how the DGP model extends to higher dimensions. Purely 
from a phenomenological perspective, it is certainly worthwhile to explore this 
wider class of infrared modifications to gravity. A second motivation, however, is 
the degravitation approach to the cosmological constant problem. The standard 
DGP model fails to exhibit degravitation. Although the DGP propagator is of 



the desired form for degravitating the cosmological constant, unfortunately (48 ) 



does not exhibit degravitation — a large brane tension generates rapid Hubble 



expansion on the brane. This is consistent with the degravitation condition ( 70 ), 
which requires a < 1/2. 

On the other hand, theories with D > 6 bulk space-time dimensions all 
correspond to a ss in the far infrared [55] . Since gravity is i'-dimensional in the 
infrared, the gravitational potential scales as l/r^~^ in this limit. Expanding 
in terms of massive states, 

$(r)=/ dMV(M^)— (71) 
Jo r r 

the spectral density must satisfy pi^'P) ^ M^^^ as AI — 0. Therefore, in the 
small momentum limit (□ — ^ 0), 

lim ^7 — ^ - / dM^— — ^ . (72) 

The integral converges for > 6. Therefore, all such theories correspond to 
a = in the infrared. For Z? = 6, the integral is logarithmically divergent, 
corresponding to rn^^^r^^) ^ log □. 
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Realizing these higher-dimensional scenarios has been notoriously difficult. 
To begin with, the simplest constructions are plagued by ghost instabilities, even 
around a flat space background |147[ I148j . Secondly, because of the higher- 
codimension nature of the brane, the AD propagator is divergent and must 
be regularized [1511 11521 1153) , usually by giving the brane a finite thickness. 
Finally, for a static bulk, the geometry for codimension N > 2 has a naked 
singularity at finite distance from the brane, for arbitrarily small tension [144] . 
(Interestingly, however, it has been argued that allowing the brane to inflate 
results in a Hubble rate on the brane which is inversely proportional to the 
brane tension for codimension N > 2 [144) .) 

Recently, it was argued that these pathologies are resolved by embedding our 
3-brane within a succession of higher-dimensional branes, each with their own 
induced gravity term, embedded in one another in a flat I?-dimensional bulk [58l 
[59l [60] . See |154j for a pedagogical review. In the simplest codimension-2 case, 
our 3-brane is embedded in a 4-brane within a flat 6D bulk, with action 



S 



Ml 
2 

2 



Mi 



bulk 



3-br 



4-br, 



d*.xV-54 (Ri + 'Cniatter) 



(73) 



The gravitational force law therefore "cascades" from AD (l/r^) to 5D 
to 61? (l/r^) etc., as we probe larger distances on the 3-brane, with the cross- 
over scales set by the ratios ms = M|/M|j and me = Mg/M|. A similar 
cascading behavior of the force law was also obtained recently in a different 
codimension-two framework |155j . Closely related work on intersecting branes 
was discussed in |156) with somewhat different motivations. See [551 11401 11571 
1158] for cosmological explorations, and |159j for self-accelerated solutions in this 
context. 

As claimed, the induced (scalar) propagator on the 3-brane is completely 
regular [58l|59]: 

Sip) - , , ^ , , (74) 



where 



+ m^gip'^) 



tanh 



(■\/""6+Py 



for p < TOg 



for p > niQ . 



(75) 



If TO5 > mg, then the Fourier transform of (74 1 for a point source yields a po 



tential that cascades from 1/r (for r <C rs) to l/r^ (for m^^ ^ r ^ "^e"^) 

(for r ^ m,g ^). Note the crucial role played by the 4-brane: in the limit 
M5 — >• 0, one recovers the logarithmic divergence, Q{p) ^ logp, characteristic 
of codimension-2 branes |151j . This can be understood intuitively as follows. 
Because of the 5D Einstein-Hilbert term, the force law must become approx- 
imately 5D at short distances on the 4-brane. Hence the lower-dimensional 
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brane behaves effectively as a codimension-l source, whose Green's function is 
therefore regular. These conclusions straightforwardly generalize to arbitrary D 
dimensions. 

The ghost issue is trickier. Perturbing around the flat space solution with 
empty branes, one finds that a ghost scalar mode propagates. In the 6D case, 
this is seen most directly from the tensor structure of the one-graviton exchange 
amplitude between conserved sources on the 3-brane in the UV limit [SH EH] '■ 



1 

□ 



T' - -Vu.T'] - -T ■ -T' . 



(76) 



We have conveniently separated terms as the sum of a massive spin-2 contribu- 
tion, with the well-known 1/3 coefficient, plus a contribution from a conformally- 
coupled scalar. The problem is with this scalar mode. The last term in (76 1 



is negative, indicating a ghost mode. This UV behavior is identical to other 
higher-dimensional scenarios [1471 1148j . 





Figure 2: Plot of the D = 7 cascading solution for the gravitational potential ^{y, z, w) 
resulting from a small tension on the codimension-3 brane. The solution is shown for ui = 
and w = 2m^^ in the case where me = m-^. In terms of the extra-dimensional coordinates y, z 
and w, the codimension-l brane is located at ui = 0, the codimension-2 brane at z = w = 0, 
and the codimension-3 brane a,t y = z = w = 0. 

However it was immediately noticed ^58] that the ghost is removed by adding 
a sufficiently large tension A on the 3-brane. (Alternatively, the ghost is also 
cured when considering a higher-dimensional Einstein-Hilbert term localized on 
the brane [1481 11491 11501 [55].) In analogy with a cosmic string in ordinary 4D 
gravity, adding tension to a codimension-2 defect leaves the induced geometry 
flat but creates a deficit angle. The exchange amplitude in this case becomes [58l 
[62] 

-4 ^ . 1 . (t;. - 1,,.T') + J T . It' . (77) 

The scalar contribution is therefore healthy provided that 

A > ^m^Mli . (78) 
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Note again the essential role played by the 4-brane kinetic term: in the limit 



M5 — >■ (keeping Mpi and Mq fixed), we have me — >■ 00, and hence (781 



cannot be satisfied. A limitation of the original derivation of (78 1 is that it was 



derived in a decoupling limit of (^73|. Recently, this result was confirmed by 



perturbing the full &D action (|73| around a background including tension on 
the 3-brane The background geometry is flat everywhere, but the extra 

dimensions show a deficit angle due to the codimension-2 source. (Just as in 



DGP, the helicity-0 and conformally-coupled scalar modes contributing to ( 77 1 
should decouple non-linearly due to the Vainshtein screening effect.) 

The D = 6 framework already exhibits degravitation: as mentioned above, a 
3-brane with tension creates a deficit angle in the bulk while remaining fiat. This 
self-tuning mechanism crucially relies on the extra dimensions having infinite 
volume: if the dimensions were compact, the brane tension would have to be 
tuned against other branes and/or bulk fluxes [160j . Since the deflcit angle 
must be less than 2tt, however, the tension allowed by the solutions considered 
in [551 EH is bounded 

A < 2t:M^ , (79) 

where the 6D Planck mass is itself constrained phenomenologically, Afg ^ meV. 
Given the geometrical origin of this bound, it is most likely an artifact of the 
codimension-2 case and is expected to be absent in higher-codimension. 

Motivated by these considerations, [ST] considered the case of -D = 7 Cas- 
cading Gravity, consisting of a 3-brane living on a 4-brane, itself embedded in 
a 5-brane, together in a 6-|-l-dimensional bulk. Including tension on the 3- 
brane, these authors derived a solution for which the induced 3-brane geometry 
is exactly flat. The hope is that in this case brane tension can still degravitate 



without being limited by a bound analogous to (79). Unlike the case of a pure 
codimension-3 DGP brane in 6-1-1 dimensions, where the static bulk geometry 
has a naked singularity for arbitrarily small tension |144j , here the bulk metric is 
non-singular everywhere (except, of course, for a delta-function in curvature at 
the 3-brane location) and asymptotically flat. This smoothness of the solution 
traces back to the cascading mechanism of regulating the propagator: the pres- 
ence of parent branes removes the power-law divergence in the AD propagator. 

Figure [2] shows the gravitational potential as a function of the 3 extra- 
dimensional coordinates, y, z and w, with the codimension-1 brane located at 
w = 0, the codimension-2 brane ai z = w = 0, and the codimension-3 brane 
at y — z ~ w — 0. Here, my = M^/Mq denotes the cross-over scale from 
6D to ID. Thanks to the cascading mechanism, which has regularized all po- 
tential divergences, the solution shown in Fig. [2] is finite everywhere. Because 
the metric depends on 3 spatial coordinates, the analytical results of [ST] was 
restricted to the weak-field approximation, corresponding to sufficiently small 
tension. Ongoing numerical work is aimed at extending these solutions to the 
non-linear regime of large tension . 
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3. Tests from Astronomical Observations 

The most stringent current tests of gravity come from laboratory and solar 
system tests, spanning millimeter to AU scales. In addition, the production of 
gravity waves predicted by GR is tested by binary pulsar observations. These 
tests are summarized above in Sec. 1 and described in Will's review [1]. 

Astrophysical probes of gravity on kpc— Mpc scales include galaxy rotation 
curves, satellite dynamics and other dynamical and lensing properties of galax- 
ies and clusters. On scales of 1 Mpc to over 1 Gpc, tests using the large-scale 
structure in the universe are being devised using the distribution of galaxies, 
dark matter (via lensing) and cross-correlations with the CMB. These astro- 
physical/cosmological tests are the subject of this section. 

Modifications in gravity can also affect the propagation of gravitational 
waves. Future gravitational wave experiments such as LISA can detect grav- 
itational waves from distant supermassive black hole pairs in the coalescence 
phase |162j ; however these tests as well as a variety of other tests of gravity in 
the strong field regime are beyond the scope of this review. 

3.1. Early Universe Tests 

The two observational probes of the early universe are Big Bang Nucle- 
osynthesis (BBN) when the universe was ^ 10 — 100 seconds old and the Cos- 
mic Microwave Background (CMB) which last scattered when the universe was 
400, 000 years old. The abundance of light elements, in particular Helium 
and Deuterium, produced during BBN is very sensitive to the expansion rate 
H of the universe and the temperature and time elapsed between the freeze out 
of neutrons and protons and BBN (see |163| for a review) . Its sensitivity to H 
and T means that the light element abundance tests the Friedmann equation, 
which can be written as: ~ Gg^,T'^, assuming radiation domination and ne- 
glecting spatial curvature. Given our knowledge of , the number of relativistic 
species at BBN, this test is usually stated as a constraint on the deviation of 
the Gravitational constant G from Newton's value Gn- Current constraints on 
G are somewhat better than 10%. The value and constraints on G are quite 
sensitive to our knowledge of (and generally assume the current values for 
non-gravitational constants such as the fine structure constant). 

The CMB at the epoch of last scattering, z ~ 1100, does not yet offer as 
clean a test of the Friedmann equation as BBN. The reason is that there is 
a significant degeneracy between G and H unless precise information on the 
CMB anisotropy power spectrum is available |164j . A higher value of G leads 
to a delay in the epoch of recombination thus causing greater damping of the 
CMB power spectrum at high angular wavenumber £. Using present data on the 
CMB temperature anisotropy power spectrum from the WMAP measurement 
plus smaller scale data |165j obtain constraints on G at nearly the 10% level. 
Combining the constraint with BBN, assuming a constant G between BBN 
and z ~ 1100, improves constraints by a factor of 3. The CMB constraint is 
limited by the fact that other parameters, in particular the spectral index of 
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the primordial power spectrum, can mimic the change in the damping due to 
deviations in G. 

With the expected measurement of the polarization power spectrum from 
the PLANCK satelhte, this degeneracy gets broken and the CMB constraints 
approach the level of 1% |165) . Thus with PLANCK and future cosmic variance 
limited measurements of the CMB polarization the constraints on G approach 
the level of current laboratory measurements. 

3.2. Large- Scale Structure Tests 

Structure formation in modified gravity in general differs [1661 11671 11681 
[T691 [TTOl [T7T1 [1721 [ml [T741 [T75l [T76l ITTTl ITTSl [MOl [T57l [158] from that in 
GR. Perturbative calculations at large scales have shown that it is promising to 
connect predictions in these theories with observations of large-scale structure 
(LSS). On small scales even larger deviations may occur around galaxy and 
cluster halos, depending on how the theory transitions to GR. Nevertheless, 
carrying out robust tests of MG in practice is challenging. The two approaches 
that have been taken are to either constrain the parameters of a particular 
model by working out its predictions for the growth of structure, or to define 
effective parameters in the spirit of the PPN formalism used to test GR in the 
solar system. While both approaches have their limitations, we shall see below 
that there has been recent progress and upcoming surveys offer great promise 
for tests of gravity. 

There are three regimes for the growth of perturbations: the superhorizon 
regime, the quasi-static Newtonian regime of linear growth, and the small scale, 
nonlinear regime. The quasi-static Newtonian regime is valid for non-relativistic 
motions and length scales sufficiently smaller than the horizon. In this regime 
(discussed in the next sub-section) the linearized fluid equations in expanding 
coordinates are sufficient to describe perturbations. In the nonlinear regime, 
while gravity is still in the weak field limit, density fiuctuations are no longer 
small and in addition the density /potential fields may couple to additional scalar 
fields introduced in MG theories. The nonlinear regime is therefore the hardest 
to describe in any general way as the nature of the coupling to scalar fields is 
theory specific. It may well be the most discriminatory for some theories as 
there can be a rich phenomenology ranging from galaxy cluster to solar system 
and laboratory scales. 

3.2.1. Metric and fluid perturbations 

In the Newtonian gauge, scalar perturbations to the metric are fully specified 
by two scalar potentials and $: 

ds^ = -(1 + 2^-) dt^ + (1 - 2$) a^{t) Ax^ , (80) 

where a{t) is the expansion scale factor. This form for the perturbed metric is 
fully general for any metric theory of gravity, aside from having excluded vector 
and tensor perturbations (see |179j and references therein for justifications). 
Note that corresponds to the Newtonian potential for the acceleration of 
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particles, and that in General Relativity $ = ^' in the absence of anisotropic 
stresses. 

A metric theory of gravity relates the two potentials above to the perturbed 
energy-momentum tensor. We introduce variables to characterize the density 
and velocity perturbations for a fluid, which we will use to describe matter and 
dark energy (we will also consider pressure and anisotropic stress below). The 
density fluctuation S is given by 

si.M)^em^, (81) 

where p{x, t) is the density and p{t) is the cosmic mean density. The second 
fluid variable is the divergence of the peculiar velocity 

0„ = v,ro7(p + p) - v-«, (82) 

where v is the (proper) peculiar velocity. Choosing 6^, instead of the vector 
V implies that we have assumed v to be irrotational. This approximation is 
sufficiently accurate in the linear regime for minimally coupled MG models. 

In principle, observations of large-scale structure can directly measure the 
four perturbed variables introduced above: the two scalar potentials 5* and $, 
and the density and velocity perturbations specified by 5 and 9^. As shown be- 
low, these variables are the key to distinguishing MG models from GR plus dark 
energy. Each has a scale and redshift dependence, so it is worth noting which 
variables and at what scale and redshift are probed by different observations. 
It is convenient to work with Fourier transforms, such as: 

5{k,t) = / d!^x 5[x,t) e"'^-^ . (83) 



When we refer to length scale A, it corresponds to a statistic such as the power 
spectrum on wavenumber k = 27r/A. We will henceforth work exclusively with 
the Fourier space quantities and drop the* symbol for convenience. 

The evolution of perturbations can be calculated in the linear regime. We 
follow the formalism and notation of jlSOj . except that we use physical time 
t instead of conformal time. We are interested in the evolution of perturba- 
tions after decoupling, so we will neglect radiation and neutrinos as sources of 
perturbations. 

The superhorizon regime is the most constrained and simplest to describe. 
As pointed out by jl79| , any metric theory of gravity that also obeys the Einstein 
equivalence principle must satisfy a universal evolution equation for metric per- 
turbations. In the conformal Newtonian gauge, for adiabatic initial conditions, 
this evolution is given by: 

^- ^^ + H<i! +\2H- = Q. (84) 

The above equation is equivalent to Eqn. (7) of |49] which uses a different time 
variable and the opposite sign for 5*. By treating the ratio of metric potential 
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as a constant parameter of a MG theory, one can solve this equation for 
a given background solution H{t). The ISW effect discussed below extends to 
very large scales and probes the superhorizon regime. 

3.2.2. Quasi-static Newtonian Regime 

In what follows, we will for the most part make the approximation of non- 
relativistic motions and restrict ourselves to sub-horizon length scales. One 
can also self-consistently neglect time derivatives of the metric potentials in 
comparison to spatial gradients. These approximations will be referred to as 
the quasi-static, Newtonian regime. Using the linearized fluid equations, the 
evolution of density (or velocity) perturbations can be described by a single 
second order differential equation: 

S + 2HS+^ = 0. (85) 

With S{k, t) ~ ^initial (fc)^(fc, t), wc Can substitute for "if in terms of 6 using the 
Poisson equation. Here we write the Poisson equation with the sum of potentials 
on the left-hand side, as this is convenient for describing lensing and the ISW 
effect. Using the generalized gravitational "constant" Gcs we have 

k^{^ + $) = -87rGcff (/c, t)p6 . (86) 

Using the two equations above, we obtain for the linear growth factor D{k,t): 

b + 2HD-jP^^-pD = Q. (87) 
(1 + $/*) ^ ' 

From the above equation one sees readily how the combination of GcS and 
<^/^ alters the linear growth factor. Further, if these parameters have a scale 
dependence, then even the linear growth factor D becomes scale dependent, a 
feature not seen in smooth dark energy models. 

We will use the power spectra of various observables to describe their scale 
dependent two point correlations. As an example, the 3-dimensional power 
spectrum of the density contrast (5(fc, z) is defined as 

{5{k, z)5{k\ z)) = (2^)3^0 (fc + k')Pss{k, z) , (88) 

where we have switched the time variable to the observable redshift z. The power 
spectra of perturbations in other quantities are defined analogously. We will 
denote the cross-spectra of two different variables with appropriate subscripts. 
For example, Ps^i denotes the cross-spectrum of the density perturbation 5 and 
the potential 

Figure [3] shows the linear and nonlinear power spectra Pssik, z) for f{R) and 
the normal branch DGP models 11811 . The dotted curves show the fractional 
departures of the f{R) linear power spectrum (left panel) and normal branch 
DGP (right panel) to ACDM. These are simply the ratio of the square of the 
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Figure 3: Power spectra for f{R) (left panel) and DGP (right panel) theories. The fractional 
deviation from A-CDM are shown for the present day linear and nonlinear power spectra 11811 . 
At high-fc (small scales), nonlinear gravitational clustering and the screening of massive halos 
alters the power spectrum. The limited resolution of the simulations precludes a definitive 
statement about the high-fc limit of the power spectra. 

linear growth factors. The strong scale dependence is evident in the f{R) case. 
In both cases significant deviations are evident at wavenumbers k }t O.l/i/Mpc. 
The DGP model shown here is the normal branch, which has stronger gravita- 
tional forces than GR (the self-accelerating branch has weaker forces than GR). 
The nonlinear power spectra have smaller deviations from GR at high-fc in part 
because the nonlinear contributions are similar for the MG models and ACDM 
at sufficiently high wavenumbers. The screening mechanism within massive ha- 
los is the other reason the deviations decrease at high-fc, but the simulations 
used do not have the resolution required to fully evaluate this contribution. 

It has been argued that a scale dependence in the MG parameters is gener- 
ically expected for scalar-tensor theories |182) . Provided the interaction can be 
expressed in the Yukawa form, one expects a correction to the growth factor 
that is quadratic in wavenumber k. The authors of [183] have also argued that 
the junction condition across the brane in DGP like theories lead to a linear 
scale dependent term that can be significant on scales approaching the horizon. 

We can also use the relations given above to obtain the linear growth factors 
for the velocity and the potentials from D. The growth factor for the veloc- 
ity divergence is given in the following sub-section, while the Poisson equation 
determines the evolution of the potentials. 

3.2.3. Tests using Lensing, ISW and Dynamical Observables 

For tests of gravity, observables that rely directly on the change in energy or 
direction of photons are distinct from those that measure the clustering or dy- 
namics of tracers such as galaxies or galaxy clusters, which have non-relativistic 
motions. We summarize the primary observables that provide tests of gravity 
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on cosmological scales in this sub-section, and consider galaxy and cluster scale 
tests in the following sub-section. We closely follows the treatment of Jain & 
Zhang [184]. 

The CMB power spectrum at angular wavenumber I is given by a projection 
along the line of sight: 



where r denotes the comoving angular diameter distance and the spherical Bessel 
function ji is the geometric term through which the CMB power spectrum de- 
pends on the distance to the last scattering surface. The function -Fcmb com- 
bines several terms describing the primordial power spectrum and the growth 
of the potential up to last scattering. We will regard Fcmb as identical to the 
GR prediction since we do not invoke MG in the early universe. (See Section 
3.1 for a discussion of CMB constraints on G and its degeneracy with other 
parameters) . 

The ISW effect: The CMB anisotropy does receive contributions at red- 
shifts below last scattering, in particular due to the integrated Sachs Wolfe 
(ISW) effect 185 . In the presence of dark energy or due to modifications in 
gravity, gravitational potentials evolve in time and produce a net change in the 
energy of CMB photons: 



The ISW effect, like gravitational lensing, depends on and probes the combi- 
nation -f $. The ISW signal is overwhelmed by the primary CMB at all 
scales except for a bump it produces at the largest scales in the CMB power 
spectrum. For this reason, it is more effectively measured indirectly, through 
cross-correlation with tracers of large scale structure at low redshift. The re- 
sulting cross-correlation signal is a projection of ^^(-4,^$) {k, x), the cross-power 

spectrum of (^'-f $) and galaxies (or other tracers of the LSS such as quasars or 
clusters). By cross-correlating the CMB temperature with galaxy over-density 
6g , the ISW effect has been detected at < 5a confidence level |186| and provides 
independent evidence for dark energy, given the prior of a spatially flat universe 
and GR. It has also provided useful constraints on MG theories as discussed 
below in Section 3.3. 

Gravitational Lensing: Lensing observables are the result of coherent 
deflections of light by mass concentrations. For the metric of Eqn. [80| the flrst 
order perturbation to a photon trajectory is given by (generalizing for example 



the gradient transverse to the unperturbed path. This gives the deflection angle 




(89) 




(90) 



Eqn. (7.72) of [T57] ): 
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Figure 4: Examples of the shear-shear and galaxy-shear power spectra for the DES (left 
panel) and a Stage-IV survey similar to LSST (right panel) |190l . The upper (green) curves 
show the galaxy-shear cross power spectrum Cgn, with foreground galaxies at 2 = 0.4 and 
background galaxies at z = 1. The lower two curves show the shear-shear power spectrum 
Ckk with two choices of redshift bins as indicated. The error bars include the sample variance 
and shape noise for the two surveys (see Table 1 for survey parameters). The shape noise 
contribution to for z = 1 is shown separately as well (dashed lines). 



formula 



(92) 



where s = qX is the path length and is the i— th component of the deflection 
angle (a two-component vector on the sky). Since all lensing observables are 
obtained by taking derivatives of the deflection angle, they necessarily depend 
only on the combination (to first order in the potentials). The convergence 

for example is given by the line-of-sight projection: 



<0) 



dz r{z)r{zs,z) „ 



H{z) r{zs) 



V^(* + $) 



(93) 



We take the sources to lie at redshift Zg . The primary observables used in weak 
lensing are the two-point correlations of the observed shapes of galaxies and the 
cross-correlation of foreground galaxies with the shapes of background galaxies. 



The metric potentials are related to the mass distribution by ( 86 1 so the 



lensing power spectra can be expressed in terms of the three-dimensional mass 
power spectrum Pss{k^z). In the small-sky-patch limit the Limber approxima- 
tion |188j gives 



dz 



H{z)a- 



■ C'(fc, z)Pss{k, z)Wl{z, Zi)WL{z, z,)lM) 



where the function C contains the modified gravity parameters: — G/Gn[{^ + 
\I/)/$] expressed in Fourier space. The lensing weight function WL{z,Zk) de- 
pends on the geometry and the redshift distribution of lensed galaxies. The 
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three-dimensional wavenumber k is given by A: = l/r{z). By dividing the galaxy 
distribution into bins in redshift |189j a number of auto and cross-spectra can be 
measured. The redshift dependence of these lensing spectra carries information 
about the growth of structure that can test gravity theories. The galaxy-shear 
cross-spectrum C^^ can be defined in a similar way to C^k^ it is proportional 
to b flmC, where h is the galaxy bias parameter. Cg^ is easier to measure and 
can be used to test gravity as discussed below. Examples of the two lensing 
spectra are shown in Figure |4] for two different survey parameters |190j . Sta- 
tistical errors are shown for the different power spectra - it is evident that if 
systematic errors can be controlled, upcoming surveys will provide percent level 
measurements |191j . 

Lensing observables probe the sum of metric potentials - this follows from 
the geodesic equation applied to photons and is therefore true for any metric 
theory of gravity. Moreover the relation of the sum of metric potentials to the 
mass distribution is very close to that of GR in scalar-tensor theories that we 
have considered [192] . Since the Einstein frame is obtained through a confor- 
mal transformation which cannot alter null geodesies, the scalar field does not 
directly alter the geodesies of light rays. Thus the deflection angle formula (92 1 
and Poisson equation ( 86 1 in the form given above are essentially unaltered in 
these scalar-tensory gravity theories. Masses of halos inferred from lensing are 
the true masses. Therefore tests of gravity that rely solely on lensing mea- 
surements can only be done by using multiple source redshifts to probe the 
growth of structure {e.g., |193| ). It is more useful to combine lensing with other 
observations of large-scale structure. 

Galaxy and cluster dynamics: MG theories do alter the force experi- 
enced by Newtonian test bodies. Thus the Newtonian potential ^ differs from 
its value in GR in both DGP and f{R) theories. It can be stronger by a factor 
of 4/3 in a certain regime in /(i?) gravity (on scales between those of chameleon 
effects and the Compton wavelength of the fa field) . This corresponds to the 
ratio <I>/^' — 1/2 with the sum -I- <I> remaining unaltered as discussed above. 
Thus for a given mass distribution, significant force enhancements can occur. 
For DGP gravity, similar force enhancements occur for the normal branch. In 
particular, the dynamical masses of halos inferred from observations relying on 
the virial theorem or hydrostatic equilibrium can be significantly larger than 
the lensing (or true) masses {e.g., |184( 1194] ). However in general the ratio of 
metric potentials is scale-dependent in MG theories, and force modifications 
can depend on halo mass and environment, so it is not straightforward to use 
dynamical and lensing masses to obtain ^/"^ . 

Constraints on the Newtonian potential ^E* on small scales are obtained using 
dynamical probes, typically involving galaxy or cluster velocity measurements. 
On sub-Mpc scales, the Virial theorem for self-gravitating systems in equilibrium 
can be used to constrain ^ in galaxy and cluster halos. Velocity tracers for 
galaxies include stars and neutral Hydrogen gas within the halos and satellite 
galaxies that orbit the outer parts of halos. For galaxy clusters the tracers are 
member galaxies and the X-ray emitting hot gas, which is also mapped using 
the Sunyaev Zel'dovich (SZ) effect. For relaxed clusters the hot gas is assumed 
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to be in hydrostatic equilibrium within the gravitational potential of the halo. 
On large scales, observables depend on the linear growth factor D(t) given 



by (851, which determines the clustering of matter and is dependent only upon 
the Newtonian potential The resulting change in the mass power spectrum 
depends on how much deviates from its GR value and the duration of time 
this deviation lasts. The results for f{R) and DGP models at 2: = are shown 
in Fig. [Sj However a direct probe of ^' at a given redshift is provided by the 
distortions of galaxy clustering in redshift space. 

Redshift surveys of galaxies provide statistical measurements of clustering 
over Mpc-Gpc scales. Redshift space distortions in the galaxy power spectrum 
arise from motions along the line-of-sight - on large scales these are sensitive to 
the linear growth factor for 9^, denoted here. It is related to D, the linear 
density growth factor via the continuity equation as: 

Dy(xaD = al3HD; P = d\nD / dlna . (95) 

Thus redshift space measurements can probes the rate of growth of structure. 
Different combinations of the information from weak lensing and redshift space 
galaxy clustering have been used to forecast tests of MG models |1951 11961 11971 
[182, 198, 190 . 

The line-of-sight component of peculiar velocities cause the observable redshift- 
space power spectrum Pgf {k, /ife) to be 'squashed' along the line of sight on large 
scales (in the linear regime) and to produce pronounced 'finger-of-God' features 
on small scales (in the nonlinear regime) [1991 1200j . The directional dependence 
of 

Pgg is given by = k\\/k, which depends on the angle between a wave vec- 
tor k and the line-of-sight direction. Although the picture is more complicated 
in reality (see j204| for a detailed discussion), it is a good approximation to 
decompose the redshift space power spectrum in terms of three isotropic power 
spectra relating the galaxy overdensity 5g = bS and peculiar velocities v: the 
galaxy power spectrum Pgg{k), the velocity power spectrum Pvv{k) and the 
cross power spectrum Pgv{k) as follows [199j 

Pg(^)(fc,Mfe) = [Pgg{k) + 2filPg,{k)+fitPUk)]F{k'filal), (96) 

where the term F {k'^ fi^a'^) describes non-linear velocity dispersion effects. The 
angular dependence in the above equation allows us to obtain the component 
power spectra from Pgg^ ■ Pgg{k), the real space power spectrum of galaxies, is 
the easiest to measure but its interpretation requires knowledge of galaxy bias. 
The pure velocity power spectrum P^y[k) has the largest error bars, while the 
cross-spectrum Pgy{k) can be estimated more easily. 

The combination of Pgy{k), which is proportional to bDy, with the galaxy- 
shear cross spectrum Cg^, is a powerful test of gravity |195j . The former probes 
the growth factor which responds to the Newtonian potential 5*, while the 
latter probes the sum of metric potentials. The ratio is independent of galaxy 
bias. Hence with appropriate redshift binning, these spectra can constrain the 
ratio of metric potentials. A recent measurement is discussed below. 
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density: max[ln(l+5)] potential: min[Y] field: min|7y/go] 




Figure 5: Simulation slices showing the density, potential and fn field for two different 
values of fnQ 12071 . Along the line of sight through a 2D slice, the maximum of the density and 
minimum values of and fn are shown using a grayscale. In the upper panels the chameleon 
mechanism is more evident - it suppresses the deviations (from GR) in the potential gradients 
in massive structures. In the lower panel the chameleon effect is much weaker due to the choice 
of a larger value of the fno parameter. These images show qualitatively the coupling between 
the Newtonian potential and the scalar field in MG gravity theories, and how small scale 
structure can vary depending on the details of the model. 
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Figure 6: Halo mass functions for f{R) (left panel) and normal branch DGP (right panel) 
models. The fractional deviation from A-CDM are shown using analytical models (shaded 
regions) and N-body simulations (symbols with error bars) |181l . See text for a discussion of 
the two definitions of halo mass. 

3.2.4- The Nonlinear Regime: Power Spectra and Halo Properties 

The nonlinear regime can be demarcated in different ways. While it is stan- 
dard in cosmology to use the density contrast being close to unity as a rule of 
thumb, for MG we demarcate it using the breakdown of the linearized equations 
for the growth of perturbations. Of particular interest is the coupling of the den- 
sity field to a scalar field such as fji for f{R) models. A number of papers have 
recently reported simulations that include such a coupling for both f{R) and 
DGP gravity [HTl HOZl Fig. [5] shows results from such a study by 

Oyaizu et al. [207] . The density, potential and scalar field fn are shown for two 
/(i?) simulations. It is evident that the qualitative difference in the fn field 
leads to a perceptible difference in the potential due to their coupled evolution. 

The two primary results from simulations are the nonlinear mass power spec- 
trum and the mass function of galaxy and cluster sized halos [1811 12091 [55] . For 
specific models these predictions allow for comparisons to data. Figures [3] and [6] 
show the measured power spectra and mass functions from N-body simulations. 
Simulations and analytical studies '194' '208', 12091 also show the transition from 
a modified gravity regime to GR inside massive halos, due to the Vainshtein or 
chameleon mechanisms for DGP and f{R) theories respectively. Thus we have 
the possibility of comparing these MG theories to observations even within and 
around galaxy and cluster halos. The advantage of such tests is that order unity 
deviations are expected. 

The mass function of cluster sized halos can show strong deviations due to 
enhanced gravitational forces from the coupling with the scalar field in MG 
theories. Fig. |6] shows the fractional deviations in the mass function for f{R) 
and DGP models (similar to the power spectra shown in Fig. [s]) |181| . For /(i?) 
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models significant departures from GR may occur if scales involved in forming 
clusters are smaller than the Compton wavelength of the scalar fn but larger 
than the scale of chameleon effects that screen halos from the modified forces 
and drive the theory to GR. Theoretical predictions require careful treatment 
of spherical collapse in MG theories as screening effects are dependent on halo 
mass and environment [194] . 

For cluster masses approching 10^^ Mq, the deviations are significant, from 
enhancements of tens of percent for f{R) to about a factor of 2 for DGP models 
(note that the DGP model shown in Fig. [6] is the normal branch DGP, which 
has enhanced forces, similar to f{R) gravity). These departures are driven by 
the deviation in the linear growth rate on a scale of ~ 10 Mpc, coupled to the 
exponential dependence of the mass function at the high mass end to this growth 
rate. In addition, the mass inferred from dynamical measurements differs from 
the true (lensing) mass for unscreened halos; this amplifies the deviations in the 
mass function as shown by the upper set of curves in Fig. [6j The resulting 
observational constraints are summarized next. 



3.3. Summary of Current Astrophysical Constraints 

Current tests of gravity on astrophysical scales rely on combining probes 
of the expansion history with the growth of structure. In the literature so far 
the observations that have been most effectively used for the distance-redshift 
relation are: CMB, SNIa and BAOs; and for the growth of structure: the CMB 
power spectrum, ISW cross-correlations, weak lensing shear correlations and the 
abundance of galaxy clusters (in the nonlinear regime). 

Perhaps the tightest formal constraints on gravity have been obtained on 
galactic scales. |2U5] and [206] used a combination of strong lensing observations 
in SDSS galaxies and the dynamics of stars to constrain the ratio of metric 
potentials. They find the ratio to be consistent with unity, to better than 10%. 

Cosmological parameter analyses have been performed by several authors to 
test gravity in the linear regime of perturbations EM EIH EEl [2131 EH 
1215112161 1217] . In the linear regime (scales >10 Mpc) current tests show consis- 
tency with GR. For the DGP scenario, |214j have obtained stringent constraints 
on the self-accelerating branch DGP model. They also find that both normal 
and self-accelerating branches of DGP require an explicit cosmological constant 
to be viable. However, generalized braneworld models such as cascading gravity 
evade these observational constraints. 

On scales above 10 Mpc a very recent test of GR [218] comes from com- 
paring galaxy-velocity and galaxy-shear cross-correlations (discussed above in 
Section 3.2.3) from the SDSS. The Eq ~ CgK./Pgv parameter introduced by 
|195j is a ratio of cross-correlations that is essentially independent of galaxy 
bias for suitable choices of galaxy redshift distributions. It is estimated to be 
~ 0.4, consistent with its value in GR: Eq = ftm{z = 0)/f3{z), where f3 is the 
logarithmic rate of growth parameter introduced above in Eqn. 95 The 20% 



level measurement of Eq by [218] spans scales of 10-50 Mpc at redshift z ~ 0.3. 

The abundance of galaxy clusters from X-ray observations has been used 
to constrain the growth factor [212^ and specific f{R) models |211j . The mass 
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Table 1: Parameters of imaging and spectroscopic surveys: sky coverage /s]jy in sq. degs., 
galaxy surface density rig"^ per sq. arcmin., the three-dimensional galaxy number density fig 
per Mpc'^, the mean redshift (z), and the maximum redshift 2max. 

function can be significantly enhanced at large masses for f{R) and DGP models 
(see discussion in the previous sub-section). Using information on the mass 
function requires nonlinear regime model predictions that include chameleon 
dynamics, hence constraints are specific to particular models. |211| constrain 
the Compton wavelength to be smaller than ^ 50 Mpc, or equivalently, for the 
present day field amplitude to be fuo < 2 x 10^^ for their f{R) model. A more 
recent analysis by [219] combines large-scale structure information with galaxy 
cluster abundances to find comparable constraints. 

Thus current tests of gravity find no indications of departures from GR. The 
tests in the literature so far are not very stringent in the linear regime, where 
model-independent constraints on MG parameters are feasible. Tests on smaller 
scales have been used to test GR by measuring the ratio of dynamical to lensing 
masses. These tests have constrained the ratio of metric potentials at the tens of 
percent level. Parameters of specific models of f{R) and DGP gravity have also 
been constrained. These tests are restricted to narrow ranges in mass/length 
scale and redshift, and are far less stringent than laboratory or solar system 
measurements. These limited tests are nevertheless a sign of progress in both 
theory and observational analysis - a decade ago virtually no tests of GR were 
available on astrophysical scales. We discuss next the prospects for tests of 
gravity from upcoming surveys. 

3.4- Upcoming Surveys 

Planned surveys that will advance low-redshift cosmological measurements 
include multi-color imaging surveys and spectroscopic surveys of selected galax- 
ies. We consider two planned imaging surveys: the Dark Energy Survey (DES) 
|220j which is expected to begin data acquisition in 2011, and a generic Stage-IV 
survey [1211 HH [Hi] whose example is the LSST survey EH]. A more 
complete list of two generations of upcoming imaging surveys, and the redshift 
range covered, is: 

• DES, PSl, HSC (2010-). < z < I. 

• LSST, JDEM, EUCLID (2016-). < z < 3. 

The surveys are characterized by sky coverage /sky, surface number density 
of lensed galaxies ni'^ and the galaxy redshift distribution. The sky coverage 
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Figure 7: Fischer constraints on modified gravity parameters from planned Stage IV surveys: 
LSST and BOSS-II [1901 . The parameters on the y-axis are qq = G/G^ and 770 = ^/^. 
The gray contours show the 68% confidence region using just imaging data while the blue 
contours combine imaging and spectroscopic (lensing and dynamical) information. The inner 
red contours also use information in the nonlinear, small-scale regime. 

for the DES is taken to be 5000 sq. degs. and for the Stage-IV survey 20000 sq. 
degs. The number density and mean redshifts of the surveys are given in Table 
[ij see |184j for the detailed redshift distribution used and other parameters. 

— • (s) • 

In order to measure the redshift-space power spectrum Pgg we consider 
spectroscopic surveys of galaxies. While primarily designed to measure the 
distance-redshift relation and H{z) using the baryon acoustic oscillations in 
galaxy power spectra, they will provide improved measurements of Pgg, Pg„, 
Pyv on large scales. Some upcoming spectroscopic galaxy surveys are: 

• LAMOST, WiggleZ, HETDEX, BOSS (2010-). < z < 3. 

• JDEM, EUCLID, BigBOSS (2016-). < z < 3. 

As an example of an upcoming survey, we consider the Baryon Oscillation 
Spectroscopic Survey (BOSS) [225^ which will target Luminous Red Galaxies 
(LRG) up to redshift z ^ 0.7 and will cover a quarter of the sky. It will obtain 
spectra of 1.5 x 10^ galaxies. In addition to the BOSS survey we consider a 
futuristic version (dubbed BOSS-II here; see also the proposed BigBOSS survey 
|226j ) with double the sky coverage compared to BOSS (to keep up with the sky 
coverage of Stage-IV survey), the same galaxy number density, and extending 
to redshift z — 1.1. Fig. [4] shows the expected accuracy in lensing spectra 
from the DES and Stage IV surveys. Fig. [7] shows the level of tests of gravity 
achievable with planned Stage IV surveys |190| using lensing and redshift space 
power spectra. The constraints on the Gravitational constant and ratio of metric 
potentials are shown using different combinations of survey data. 

There is another kind of redshift survey based on 21 cm radiation from 
neutral Hydrogen at high redshift. The most ambitious planned survey is the 
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Figure 8: Tests of gravity at different lengtfi scales. Red lines shows observations that probe 
the sum of metric potentials via weak and strong gravitational lensing or the ISW effect. Blue 
lines show dynamical measurements that rely on the motions of stars or galaxies or other 
non-relativistic tracers. With upcoming surveys, scales of 100s of Mpc will become accessible, 
while the accuracy of measurements will improve on all scales. 
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square kilometer array (SKA) |227 'f2015— ). SKA has the potential to detect 
~ 1 billion galaxies over < z ^ 1.5, with a deeper survey extending to z ~ 5, 
through 21cni line emission of neutral hydrogen in galaxies. If successful, it will 
provide high precision measurements of the distance-redshift relation through 
BAO's |228| . and tests of MG |195| through: weak lensing, velocity measure- 
ments through redshift distortions of galaxy clustering, and ISW measurements 
through CMB-galaxy cross-correlations. 

X-ray and SZ cluster surveys (e.g. ACT [229] and SPT [230]) aim to measure 
the mass function of galaxy clusters out to z ~ 1. Cosmological applications 
will depend on supplementary optical data to get photometric redshifts of the 
detected clusters. The comparison with lensing masses is essential for the mass 
calibration, but it will also provide a test of gravity by comparing dynamical 
and lensing masses. 

CMB temperature and polarization maps provide high-z constraints and also 
measurements of the ISW effect and CMB lensing, which are probes of \1/ -I- $ 
at lower redshift. In addition to the all-sky PLANCK mission from space, 
a number of surveys from the ground will provide polarization maps at high 
angular resolution. 

We have indicated the approximate redshift range over which these surveys 
will provide accurate measurements. Different observables that overlap in red- 
shift and length scale in the range z ~ 0.3 — 1 and at scale A ~ 10 to several 100 
Mpc can be compared in the linear and quasilinear regimes of structure forma- 
tion. It is likely that MG effects are significant on these scales. As emphasized 
above, on smaller scales of order a Mpc or less, even larger deviations may be 
expected. Fig. |8] shows the length scales expected to be probed at low redshift 
by lensing and dynamical probes. 

We can expect that theoretical developments will continue to influence ob- 
servational approaches since no complete formalism exists to provide model 
independent tests of MG theories on astrophysical scales. It is only for linear 
fluctuations in the quasi-static, Newtonian regime that the effective Gravita- 
tional constant and the ratio of metri potentials are useful effective parameters. 
New tests in the nonlinear regime will require detailed models and simulations 
that incorporate the screening mechanisms for specific theories. Combinations 
of laboratory and solar system tests with astrophysical tests are also promising 
probes of screening mechanisms. Finally, we note that we have not explored the 
feasibility of distinguishing modifications of gravity from general dark energy 
scenarios. Indeed, it is likely that dark energy with small scale clustering that 
allows for any desired anisotropic stress and pressure will not in practice be 
distinguishable from the kinds of MG theories we have discussed (e.g. j51 ll84p . 

4. Discussion 

We have reviewed modified gravity theories with interesting cosmological 
consequences. We focused on higher dimensional approaches and scalar-tensor 
theories, such as chameleon// (i?) and symmetron theories. We classified these 
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models in terms of the screening mechanisms that enable such theories to ap- 
proach general relativity on small scales. We described general features of the 
modified Friedman equation in such theories and how they may be distinguished 
experimentally. The second half of this review discussed experimental test of 
gravity in light of the new theoretical approaches. We introduced a variety of 
astrophysical tests that are likely to be valuable in testing gravity theories. Re- 
cent progress in modeling the growth of perturbations in MG theories and new 
observational analyses have led to the first astrophysical tests of gravity. We 
summarized current constraints, upcoming surveys and prospects for new tests. 
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